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ABSTRACT 



We present results from the X-ray portion of a multi-wavelength study of local 
ULIRGs and QSOs called QUEST (Quasar-ULIRG Evolution STudy). The data 
consist of new and archival X-ray data on 40 ULIRGs and 26 PG QSOs taken 
with Chandra and XMM-Newton. A combination of traditional and hardness ra- 
tio spectral fitting methods is used to characterize the X-ray properties of these 
objects. The absorption-corrected 2-10 keV to bolometric luminosity ratios of the 
ULIRGs and PG QSOs suggest that the likelihood for dominant nuclear activ- 
ity increases along the merger sequence from "cool" ULIRGs, "warm" ULIRGs, 
infrared-bright QSOs, and infrared-faint QSOs. The starburst dominates the to- 
tal power in ULIRGs prior to the merger, and this is followed by rapid black hole 
growth during and after coalescence. These results are in general agreement with 
those obtained in the mid-infrared with Spitzer and recent numerical simulations. 

Subject headings: galaxies: active — galaxies: starburst — X-rays: galaxies 



1. Introduction 



Since their wide-spread detection by the IRAS satellite over 20 years ago, the energy 
source in luminous (L IR > 10 11 L Q ) and ultraluminous (L IR > 10 12 L Q ) infrared galaxies 
(U/LIRGs) has been under debate. It is thought that reprocessed starlight and/or ac- 
tive galactic nuclei (AGNs) are responsible for the enormous luminosities these objects ex- 
hibit in the in f rared , but the exact contribution from each component is hard to determine. 



Sanders et al.l ( 119881 ) proposed that U/LIRGs and quasars are linked through galaxy mergers 



and ULIRGs are simply the dust-enshrouded phase in the merging process. The authors sug- 
gested that the energy source evolves along the merger sequence from starburst-dominated 
LIRGs to AGN-dominated ULIRGs and ultimately quasars. The question of the energy 
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source in U/LIRGs has cosmological implications. At higher redshifts (z > 1), U/LIRGs 
are a significant populat i on, contributing i n large part to the cosmic star formation (e.g., 
Sanders fc Mirabel! Il996t blain et al.l l2002h . About a third of the high-z ULIRGs are on- 



going gas- rich mergers (e.g., iDaddi et al.l 120071 ; IShapiro et al.ll2008l ; iForster Schreiber et al. 



200 9J). The fraction of high- z ULIRGs whi ch are m e rgers increases with luminosity (e.g. 
Tacconi et al.l 120081 ) . Merger simulations by iLi et al.l (120071 ) suggest that vigorous star for 
mation is needed to reproduce the observed dust properties of high redshift quasars. ILi et al. 
( 120081 ) were able to reproduce the observed spectral energy distribution of SDSS J1148+5251, 
a quasar at z ~ 6.4 using these merger simulations. Their work support the starburst- 
to-quasar evolutionary scenario for some luminous quasars at high redshifts. However, at 
high redshifts, secular processes may play a more dominant role than in the local universe 
and the merger scenari o may not be necessary for the formation of some quasars (e.g., 
Dekel &: Birnboimll2006l ). If the generalized merger scenario holds true, then a large fraction 
of the high- 2 ULIRGs may be the progenitors of some of the present-day quasars. 

Higher redshift objects are difficult to study, so as a first step in understanding the 
role U/LIRGs play in the formation of some quasars and elliptical galaxies, we utilize local 
U/LIRGs and quasars as laboratories for studying merger physics. Past ground-based optical 
and near-infrared photometry have shown that essent ially all ULIRGs show signs of interac - 



tion, indicating that they are ongoing mergers (e.g., iKim et al.ll2002l ; IVeilleux et al.l 120021 ). 



Recent simulations s upport the scenario that quasars can indeed be formed through gas-rich 
galaxy mergers (e.g., iHopkins et al.ll2005l 120061 120081 ). Therefore, the two populations may 
be linked. 

To further study this relationship, we are conducting a comprehensive, mult i- wavelength 
imaging and spectroscopic survey of local ULIRG and QSO mergers called QUEST - 
Quasar/ULIRG Evolution STudy. The goal of the project is to investigate the evolutionary 
link between the two groups of luminous objects in the nearby universe. The sample consists 
of z < 0.3 1-Jy ULIRGs and Palomar Green (PG) quasars. The 1-Jy sample of ULIRGs is 



a group of very well studied objects in the infrared and optical (e.g., IVeilleux et al.l 12002 



Kim et al.l 120021 ) . We pair the ULIRGs with another group of objects with a wealth of 



ground- and space-based data — the PG QSOs. At luminosities above 10 L , ULIRGs 
and PG QSOs are the only two types of extragal actic sources with simila r space densities 
and bolometric luminosities in the local universe (jSanders fc Mirabellll996|). While i n some 
respect the PG QSOs are not representative of all quasars as a group, IJester et al.l (120051 ) 
have shown that for PG QSOs with redshifts belo w 0.5, the selection bi ases are minimal. 
The detailed QUEST sample selection is described in lVeilleux et al.l (J2009a|) (hereafter V09a) 
and references therein. 
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To date, optical and infrared portions of QUEST have found that ULIRGs and PG QSOs 
have similar gala ctic structure and stellar ho st dynamics. Results from an HST NICMOS 
imaging study by IVeilleux et al.l ( 120061 . l2009bl ) suggest that AGN-like objects, including the 
QSOs, are generally of early morphological type and have less pronounced merger-induced 
morphological anomalies than systems with LINER-like or HII region-like spectral types and 
cooler infrared colors (/25//60 — 0.2). Infrared-bright QSOs generally have more pronounced 
merger-induced morphological anomalies than infrared-faint QSOs. The velocity dispersion 
distributions in ULIRGs resemble those in intermediate mass ellipicals/lenticulars with mod- 
erate rotation and the black hole masses of these ULIR Gs are estimated to be of the order 
10 7 — 1O 8 M , as demonstra ted by iDasvra et al.l (l2006al Jbr) using VLT spectroscopic data and 



the Mbh ~ <y relationship of iTremaine et al.l (120021 ). The black hole mas ses derived from sim - 
ilar data on a dozen PG QSOs agree with those of coalesced ULIRGs (IDasyra et al.l 120071 ) . 
suggesting that the bulk of the black hole growth takes place in the ULIRG phase of the 
merger. 

Ultimately, if ULIRGs and PG QSOs are linked through mergers, then we should see 
an evolution of the energy production mechanism along the merger sequence. As the merger 
ages, the dominant source of infrared radiation should change from starburst (triggered by the 
interaction) to nuclear activity. As part of QUEST, we performed Spitzer IRS observations 
of PG QSOs which showed that starbursts are responsible for at least ~30%, b ut likely 
most, of the far-infrared (FIR) luminosity of the quasars (jSchweitzer et al.l 120061 ). V09a 
found that, on average, AGN contribute ~ 40% of the bolometric luminosity in the QUEST 
ULIRGs. The AGN contribution ranges from ~15% - 35% among cool optically classified 
H Il-like and LINER ULIRGs, and ~50 and ~75% among warm Seyfert 2 and Seyfert 1 
ULIRGs, respectively. This number exceeds ~80% in PG QSOs. Thus, a trend of increasing 
AGN contribution is seen along the merger sequence, in general agreement with the standard 
ULIRG-QSO evolution scenario. 

As ULIRGs are the results of the merging of gas-rich disk galaxies, obscuration may 
be an issue for optical and infrared observations. While radio emission is less affected by 
intervening material, the bolometric luminosity in the radio band is insignificant to prove that 
accretion onto supermassive black holes is the dominant energetic process. X-ray emission, 
on the other hand, generally contributes more significantly than the radio emission to the 
bolometric luminosity of an AGN and is also less affected by obscuration than the optical 
and infrared emission. X-ray observations therefore have the potential to isolate the AGN 
contribution from that of the starburst, assuming the former arises from a compact region 
associated with the accretion disk. Chandra and XMM-Newton are excellent complementary 
instruments for this purpose. The high spatial resolution of Chandra is ideally suited to 
search for unresolved hard X-ray (2-10 keV) nuclei indicative of AGNs. However, the absence 
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of such a nucleus does not necessarily imply a starburst origin since a large gas column 
(Nh > 10 24 cm -2 ) in the line of sight (possibly due to a torus) could strongly attenuate the 
hard X-rays. In this case, X-rays emitted along the polar axis may be electron scattered into 
the line of sight, with the signature of a n Fe Kq line of large equivalent width (~1 keV; e.g., 
Ghisellini et al.lll994t iKrolik et al.lll994j ). This is where the excellent sensitivity and spectral 
resolution of XMM-Newton near the Fe K complex become most useful for AGN diagnostics. 

In the present paper, we focus on the X-ray properties of the QUES T ULIRGs and 
PG QSOs. There have been previous X-r ay snapshot studies of U/LIRGs (e.g.. lFranceschini et al 
20031 ; iPtak et al.ll2003l ; iTeng et al.ll2005l ). but the sample selection criteria for these surveys 
differed, making systematic comparisons difficult. The analyses of the X-ray data presented 
in this paper are performed in an uniform manner with an unprecedented large number of 
ULIRGs and PG QSOs. It is hoped that the contributions of the starburst and the AGN 
to the bolometric luminosity in these objects can be quantified systematically and with sta- 
tistical significance. Our current sample consists of 40 ULIRGs and 26 PG QSOs from the 
QUEST sample that have publicly available X-ray data from Chandra and/or XMM-Newton. 
These represent 51% of the ULIRGs and 79% of the PG QSOs in the full QUEST sample. 
Table ED lists the QUEST objects with X-ray data. 3C 273 and PG 0157+001 are both 
ULIRGs that are also PG QSOs. For the purpose of this paper, and to be consistent with 
the previous papers in this series (e.g., V09a), we consider these to be ULIRGs. To improve 
the statistics of our analysis and extend the range in infrared luminosity, we also add to the 
sample 26 non-QUEST U/LIRGs from the Chandra archiv e. These objects are part of the 
Revised Bright Galaxy Survey (RBGS; ISanders et al.ll2003f ) and are discussed in more detail 
in Iwasawa et al. (2010, submitted). 

The organization of this paper is as follows. The techniques used for taking and reducing 
the data are discussed in SJ2J The spectral fitting analyses for both ULIRGs and PG QSOs 
are described in JJJ We address the origin of the soft excess seen in PG QSOs in §U In §3 
we comment on the possible evolutionary link between ULIRGs and quasars, combining the 
X-ray data with recent Spitzer and HST observations. The results are summarized in §6j 
In Appendi x A, w e examine the reliability of the hardness ratio (HR) method developed in 
Teng et al.l ( 120051 ) for low-count X-ray sources. Throughout this paper, we adopt H = 70 
km s" 1 Mpc" 1 , VL M = 0.3, and fi A = 0.7. 



2. Observations and Data Reduction 



The observations of ULIRGs and PG QSOs presented in this paper are selected from the 
Chandra and XMM-Newton archives as well as our own (PI: Veilleux) guest observer (GO) 
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programs (Chandra cycle 10 and XMM-Newton cycle 7). Only ACIS-S data are considered 
for the Chandra analysis. Similarly, only EPIC data are considered for the XMM-Newton 
observations. Table [2] lists details on the available observations. 

2.1. Data Calibration and Extraction 

2.1.1. Chandra Observations 

The reduction of the archived Chandra data was performed using CIAO version 4.1.1 
and CALDB version 4.1. The cycle 10 GO data were reduced using CIAO version 4.1.2 and 
CALDB version 4.1.3. The Science Analysis Threads for ACIS dataQ outline the procedure 
used to process and reduce the data. Table [2] lists the total exposure for each observation 
after the selection of good time intervals where the data are not affected by background 
flares. 

The source spectra were extracted with circular regions centered on the source. For 
most observations, background spectra were selected from an annular source-free region that 
surrounds the nuclear extraction area. However, in particularly crowded fields, a nearby 
circular, source-free region was used. The sizes of the regions vary depending on the angular 
extent of the sources, especially in U/LIRGs where there are a lot of extended diffuse emis- 
sion. The region sizes range between 15 to 50 arcseconds and were maximized to include all 
emission from the galaxies, unless limited by nearby sources or gaps in the CCD detectors. 

2.1.2. XMM-Newton Observations 

The archived XMM-Newton data were processed using the XMM-Newton Science Anal- 
ysis System (SAS), version 7.1.0, released on 2007 July 8. The event lists were re-calibrated 
with the latest available calibration files as of 2008 July. The Cycle 7 GO data on PG 0838+770, 
PG 1435—067, and B2 2201+31A were reduced using SAS version 8.0.1. The standard pro- 
cessing procedures outlined in § 4.13 of the XMM-Newton SAS User's Guide (Issue 5.0) were 
followed for both archived and GO data. Times of high background flares were flagged and 
the total good time interval for each observation is listed in Table |2j The standard method of 
background screening involves discarding time intervals affected by background flares where 
the background count rates at energies above 10 keV are above the recommended thresholds 
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of 0.35 counts s 1 and 1 counts s 1 for EPIC-MOS and E PIC-pn data, respectively. While 



the MaxSNR method introduced by lPiconcelli et al.l (120041 ) would maximize the total net ex- 



posure times of the data, this method is only appropriate to use for data of high-flux sources 



X-ray observations of U/LIRGs (e.g., iFranceschini et al.ll2003l ; iPtak et al.l 120031 : iTeng et al. 



20051 ) have shown that these data have low signal-to-noise ratios (SNRs). Therefore, we have 
conservatively chosen to use the standard method of background screening for a consistent 
treatment of ULIRG and PG QSO data. 

The cross-calibration between the EPIC-MOS and the EPIC-pn cameras has some time 
and energy dependencies^]. Complications to the fits may result if spectra from both detectors 
are modeled simultaneously. Therefore, we chose to model the spectra from only the EPIC- 
pn camera due to its high quantum efficiency and the inability to extract background spectra 
when the small window mode was used in many of the EPIC-MOS observations. The EPIC- 
MOS data were used only in the few cases where the EPIC-pn data were unavailable or have 
much lower SNRs than the EPIC-MOS data: PG 0838+770 has very low SNRs, PG 1244+026 
has very few detected counts above ~5 keV, PG 1613+658 and PG 1626+554 have data 
highly affected by background flares. 

The SAS task EPATPLOT was used to determine whether the observations on the PG 
QSOs were affected by pile up. Only three sets of observations were found to be piled-up: 
PG 0844+349, the 17 December 2004 observation of PG 1116+215, and PG 1426+015. In 
these cases, the observations were re-extracted using annular regions to exclude the central 
part of the source which is the most susceptible to pile-up, as recommended by the SAS 
User's Guide. This pile-up correction method is reliable and does not affect the shape of the 
output source spectrum|f|. 

The XMM-Newton source spectra were also extracted with circular regions centered on 
the source ranging between 25 to 50 arcseconds depending on the angular sizes of the objects. 
For EPIC-pn data, the pointing center falls near a CCD gap. Thus, the background was 
selected from a nearby circular region in which no obvious background source resides. The 
background regions for EPIC-MOS data were selected in a similar fashion. 

Where there are multiple observations of the same source with the same instrument and 
filters, an average spectrum is created using the FTOOLS task MATHPHA giving each input 
spectrum equal weight. The sources where this was performed are the EPIC-pn observations 
of PG 0050+124, PG 1116+215, PG 1440+356, and PG 1501+106. 



2 XMM-Newton Calibration Documentation: 

3 XMM-Newton Calibration Documentation: 



http://xmm2.esac.esa.int/docs/documents/CAL-TN-0052-5-0.ps.gz 



http://xmm2.esac.esa.int/docs/documents/CAL-TN-0036-l-0.ps.gz 
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3. Spectral Analysis 

The spectral analysis was performed using XSPEC vl2.5.0. All quoted errors are 90% 
limits on one parameter (A% 2 or Ac-stat = 2.706). The errors of the derived values in the 
rest of this paper are assumed to be at the 90% confidence level. Because of the differences in 
calibration, the effective energy range is 0.3-10.0 keV for EPIC-pn, 0.6-10.0 keV for EPIC- 
MOS, and 0.5-8.0 keV for ACIS. For consistency between the different detectors and with 
the literature, the soft X-ray band measurements are made between 0.5-2.0 keV while the 
hard X-ray band measurements are made between 2.0-10.0 keV where the best-fit models are 
used to extrapolate the 0.5-0.6 keV and 8.0-10.0 keV measurements from the EPIC-MOS 
and ACIS observations, respectively. 



3.1. Results: PG QSOs 



Depending on the number of detected counts, the extracted spectra of the PG QSOs 
are binned differently for spectral modeling. Their source spectra were binned to at least 50 
counts per bin with the exception of those from fainter sources with relatively short inte- 
gration times (PG 0838+770, PG 1001+054, PG 1004+130, PG 1126-041, PG 1244+026, 
PG 1309+355, PG 1411+442, PG 1426+015, PG 1435-067, PG 1613+658, PG 1626+554, 
and PG 2214+139) which were binned to at least 15 counts per bin so that y 2 statistics 
would be applicable. 



In modeling the spectra, we take the same basic approach as in ITeng et al.l ( 120091 ): 
first, we assume a simple power-law distribution absorbed only by the Galactic column to 
describe the emission from the AGN. If the model is not a satisfactory fit to the data, 
then we consider adding a MEKAL component to describe the starburst, absorption by 
intervening material near the central source, and emission lines to model Fe Ka and lines 
of other elements, if applicable. The F-test is used to determine whether the additional 
components to the basic model are significant^. We assume P F-test < 0.001 for significant 
additions. Hereafter, we refer to these models as the Power-Law (PL) models. Table [3] lists 
the best-fit parameters of the PL models to the PG QSOs and these are shown with the 
spectra in Fi gure [D More complex models involving reflection and scattering such as those 
presented in iPiconcelli et al.l (12005) (her eafter P05 models), the blurred ionized reflection 
model presented in lCrummy et al.l (120061 ) (C06), and a constrained version of the C06 model 



4 Under certain c onditions, such as test ing; for a spectral line, it is inappropriate to use the F-test for 
model selection. See IProtassov et al.l (|2002i ) for caveats and details. 
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with fewer free parameters (C06con, described below) were also considered. 

The P05 model consists of a power law for the AGN plus four different continuum models 
to explain the soft excess often seen in AGN spectra (see § H|). These include a blackbody, 
multi-color blackbody, bremsstrahlung, and a second power law. Some of the sources also 
required additional absorption edge features between 0.6-0.8 keV. The best-fitting model 
varies for each source and according to their study, there does not seem to be an universal 
model for the quasar spectra. On the other hand, C06 selected a relativistically blurred 
reflection model as their universal model for the quasar continua. The model is a relativistic 
convolution (KDBLUR in XSPEC) of a photoionized disk reflection model (REFLION). 
The model assumes that a semi-infinite slab of optically thick cold gas of constant density 
is illuminated by emission from the cor ona. The r eflected and the direct components are 



then convolved with a Laor line profile (ILaorl Il99ll ) to create the effect of blurring from a 



relativistic accretion disk. The blurred reflection model was invoked to explain the smooth 
soft excess features and the lack of strong iron lines seen in the PG QSOs. Unlike the P05 
models, the C06 model is an universal model and can thus be applied to all the quasars in our 
sample. We first checked that we have applied the C06 model properly by ensuring that the 
best-fit values are consistent with those published in C06 for observations that overlap both 
samples. We then applied the C06 model more widely to all of our quasar spectra. The range 
of best-fit parameters derived from the C06 models require extreme values for the accretion 
disk. For example, the disk emissivity index (e, the power law dependence of the emissivity, 
r~ e ) ranges from 1.3 to 10.0, the upper limit of allowable values. Therefore, we constrained 
the C06 model by fixing some of these parameters at their more widely accepted values to see 
if it would yield acceptable fits. In our constrained version (C06con), the emissivity index of 
the disk in the KDBLUR model component is fixed at the more commonly accepted value 
of 3.0. Since the iron abundance is less reliable, we conservatively fix the abundance of the 
REFLION component at solar. For the majority of objects, the inclination of the disk is 
also fixed at the default value of 30 degrees in order for the fit to converge so that xl < 2.0. 
We find that the C06 models perform significantly better than the C06con models in only a 
few instances. 

TableslUEJ and El list the best-fit parameters to the P05, C06, and C06con models to the 
PG QSO data, respectively. Figure [2] shows comparisons of the reduced x 2 values between 
the PL and the P05 and C06 models; the PL models appear to be as good as, or better than, 
the P05 and C06 models in fitting the data. In terms of the basic spectral properties, the 
PL models seem to be equally good in determining the photon index as both P05 and COG. 
Figure [3] is a comparison of the photon index as determined by the three methods. The P05 
and PL models appear to be consistent with each other; however, the PL model requires 
softer spectra than the C06 models. Nevertheless, the photon indices as determined by all 
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three methods are within the range previously measured in the PG QSOs by other authors 
(e.g., 



Porquet et al.ll2004f > 



As the figures show, with the exception of a handful of observations, none of these 
models is statistically favored. However, variability in quasar spectra is well known, so 
the ideal model of the data should also naturally explain the observed variability seen in 
these objects. Thus, based on variability arguments, the C06 model appears to be the least 
likely explanation for the quasar spectra. C06 evaluated only single-epoch quasar spectra 
from XMM-Newton. They chose the longest available observation for objects with multiple 
data sets in the archive. The PL models can easily explain the variability by a change in 
the absorbing column (see Table [3]); physically, t his can be due to a clumpy torus viewed 
approximately edge on (see also iTeng et al.l 120091 ) . In order for the C06 model to explain 
the variability (particularly the significant variability of PG 0050+124 and PG 1501+106), 
large changes in T, the emissivity index, and the disk inclination are required (see Table [5]). 
While it can be argued that the changes in T and emissivity index are due to changes in the 
accretion flow of the black hole, the required change of about 30° in the observed inclination 
of the accretion disk over a the period of a few years is hard to justify. 



3.2. Results: ULIRGs 



As with previous U/LIRG surveys (e.g.. lFranceschini et aU2003l ; |Ptak et al.ll2003t ITeng et al 



20051 ). most ULIRGs in the present study were detected, but with many fewer counts than 
the average PG QSO. We separated them into three "brightness" categories: weak, moderate, 
and strong for the purpose of the analysis. 

The "weak" sources are those with count rates < 0.01 and 0.05 counts per second 
when observed by Chandra and XMM-Newton, respectively. With the time allocated, these 
sources do not have enough counts for traditional spectral fitting and th eir spectral proper ties 



and fluxes are measured using the HR method (see Appendix A and ITeng et al. 



2005, for 



details). The spectral properties derived using the HR method for the weak ULIRGs are 
listed in Table [7J 

The "moderate" sources are those with relatively low count rates, but the exposure 
times are long enough to obtain low signal-to-noise spectra with more than 100, but less 
than 1000, counts. For these sources, their spectral properties are modeled using both the 
HR method and traditional spec tral fitting . The traditional spectral fits were performed 
using the Cash Statistics (c-stat ICashl 1 19781 ) option in X SPEC on u n binne d spectra. The 
details of the fitting procedure are the same as those in ITeng et al.l (120081 ) . Only the PL 
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models are applied to these low signal-to-noise spectra. The best-fit properties are presented 
in Table EJ 

Finally, the "strong" sources are those with relatively high count rates and relatively 
high signal-to-noise spectra having more than 1000 counts. The spectral modeling was 
performed using the \ 2 statistics option in XSPEC as the spectra were binned to at least 
15 counts per bin, with the exception of 3C 273, which was binned to at least 100 counts 
per bin. Again, only the PL models were applied to these ULIRG spectra. The spectral 
properties of these sources are also listed in Table 

In general, the 2-10 keV band is less likely to be affected by obscuration and this is 
where AGN emission, if present, would dominate. Figure H] is a comparison of the 2-10 keV 
and infrared luminosities of the U/LIRGs and PG QSOs. The U/LIRGs include the 26 
RBGS objects from the Chandra archive for a total of 66 objects. For AGN-dominated 
objects like the PG QSOs, the 2-10 keV luminosity is consistently about 3% of the infrared 
luminosity. The 2-10 keV luminosity of the U/LIRGs, however, vary between 0.0001 to 3% of 
the infrared luminosity with the majority of the ULIRGs falling between log(L 2 -iofceV / Lir) 
of (-4.5, -1.5). Comparison with the effective mid-infrared optical depth, r e ff, derived by 
V09a suggest that AGN-dominated objects also have the smallest extinction. 

We also compared the 2-10 keV luminosity with the bolometric luminosity of the 
U/LIRGs and PG QSOs; the 2-10 keV to bolometric luminosity ratio is our proxy for 
AGN dominance. The bolometric luminosity is defined as L& / = llhL IR for the U/LIRGs 
and Lboi = 7L ° + L IR for the PG QSOs. The left panel in Figure [5] plots this ratio 

5100/1 

against the total bolometric luminosity of our sample. For the U/LIRGs there is simply a 
small shift in both axes from Figure HJ The right-hand panel shows the same ratio plotted 
against the AGN contribution to the bolometric luminosity as derived by V09a. Again, we 
observe a spread of three orders of magnitude in the 2-10 keV to bolometric luminosity for 
the U/LIRGs. The importance of nuclear activity relative to that of the starburst increases 
slightly with the bolometric (infrared) luminosity of the ULIRG. We find again that the less 
obscured objects are more likely to be AGN-dominated. 



Veilleux et al.l (l2009bl ) derived black hole mass estimates us i ng th e if-band elliptical host 



magnitude to black hole mass relation from iMarconi fc Huntl (120031 ) . The derivation does 



not include dust extinction beyond the nuclear regions of the hosts; this may possibly lead 
to underestimates of the black hole masses. The presence of recent /ongoing non-nuclear star 
formation is also not excluded which may cause overestimates in the black hole mas s es. Ap- 



plying these photometrically derived black hole mass estimates from I Veilleux et al.l (J2009bJ) 
to our X-ray observations, we calculate the 2-10 keV-to-Eddington luminosity ratio, our 
proxy for the Eddington ratio. Figure [6] compares the Eddington ratio derived from the 
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X-ray methods with that from the mid-infrared methods presented in V09a. Since the X-ray 
values are only the 2-10 keV luminosity and not the bolometric luminosity of the AGN, 
there is a shift of a factor of about 30 - 100 between the horizontal and vertical values. The 
two methods are linearly correlated for the PG QSOs, but a discrepancy between the X-ray 
and mid-infrared AGN diagnostic methods is observed for the U/LIRGs. The cause of this 
discrepancy is discussed in § 15.21 



4. An Universal Spectral Model for the PG QSOs 



PG QSOs have been very well studied in X-rays and at other wavelengths. However, 
many X-ray studies on PG QSOs in the literature have focused on only a handful of ob- 
jects. Most studies have found that the X-ray spectra of PG QSOs tend to be more or less 
featureless, with sometimes small emission features near the iron K complex arising from 
neutral or ionized iron. Often, these iron emission lines have relatively small equivalent 
widths (~100 eV) and are generally narrow. The continua of the X-ray spectra can be fit by 
a power law with T ~1.8 for the AGN, but the X-ray spectra of several PG QSOs (as well as 
many Seyfert galaxies) also contain what is termed the "s o ft exc ess" where there is emission 
in excess of the power-law below ~2 keV. iBianchi et al.l (120081 ) uses a model-independent 
method of comparing the 0.5-2 keV and 2-10 keV flux ratios as an independent measure of 
the soft excess strength. The authors confirm that the soft excess is a common feature in 
low-redshift active galaxies. 

The presence o f the soft excess in active galaxies was first identified in the 1980s (e.g., 
Arnaud et al.1119851 ). Ever since, the origin of th e emis sion has remained a puzzle. In recent 
years, several authors including iPorquet et al.l ( 120041 ) . P05, C06 have performed system- 
atic analyses on large groups of quasar spectra in an attempt to find the best model that 
describes the exce ss emission. The cur rent leading models for the soft excess include Comp- 



ton up-scattering (IPorquet et al. 



( IGierlinski &: Don 



2004 



2004 ), blurred refl ection model (C06), complex absorption 



Sobolewska fc Pond 120071 ). and other models such as discussed in 



P05. 



Porquet et al.l (120041 ) at first thought that the soft excess originates from the inner 



accretion disk. This model explains naturally the smooth transition from UV accretion disk 
emission to the soft excess. However, the inferred blackbody temperatures from the modeling 
are too hot to be explained by direct emission from a thin accretion disk for reasonable values 
of the black hole mass. Thus, they prefer their alternate model, the Compton up-scattering 
of the extreme ultraviolet photons from the accretion disk to form the soft excess. 
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C06 presented a blurred reflection model as a universal model for the PG QSO spectra. 
Blurred reflection, caused by the relativistic motion in the accretion disk, is invoked due to 
the lack of broad iron lines observed in these objects. While this model fits the spectra well, 
it also requires extreme values for some of the model parameters, as explained in § 13.11 



Gierlinski fc Pond ( 12004 ) suggested that the soft excess is caused by a broad absorp- 



tion trough at ~2-5 keV, possibly related to an accretion disk wind. The smoothness of 
the observed soft e xc ess, however, canno t be r eproduced in the latest simulations (e.g., 
Schurch et al.l 120091 ) . ISobolewska fc Pond (120071 ) favor a complex absorption model based 
on their spectral modeling of two AGNs (PG 1211+143 and 1H0707-495) with large observed 
soft excesses. 



Instead of a universal model for the soft excess, iPiconcelli et al.l (120051 ) found in their 
survey of 42 PG QSOs that the X-ray continua of these quasars are well fit by a com- 
bination of four different models. The four models are blackbody, multicolor blackbody, 
bremsstrahlung, and power law. 

From a statistical stand point, all of the above models (as well as the PL models pre- 
sented in § !3.ip are equivalent. The fitting statistics of the models to the data are very similar. 
Below, we describe a new approach to help us discriminate between the various models for 
the soft excess, where we model all of the PG QSOs with a single universal model, varying 
only a few key parameters from object to object. We do not know the actual origin of the 
soft excess, so instead of making assumptions about which model describes physics that we 
do not fully understand, we look at the data in a more general way. 

The soft excess is a common feature in low-redshift active galaxies, including quasars. 
Therefore, it can be assumed that all QSO spectra have the same basic shape arising from 
the same physical phenomena. Let us further assume that there is a single (yet unknown) 
ideal model that describes the quasar spectra in the form of a power law (for the AGN) plus 
an additional component. Thus, each observed quasar spectrum is then a random variation 
of the basic ideal model. If we fit all the spectra simultaneously with the same model, then 
we should be able to define a median model for the quasars as a class. This median model, 
then, should help us identify a favored model that best fits all objects. 



4.1. Modeling PG QSOs as a Class 

For this multi-source fitting, we chose to use only XMM-Newton EPIC-pn observations 
to minimize cross-calibration issues. This includes 22 objects except PG 0804+761 (unde- 
tected), PG 1244+026 (few counts above ~5 keV), PG 1613+658, and PG1626+554 (no 
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pn data). Average spectra for objects with multiple observations were created using the 
FTOOLS task MATHPHA and used for the multi-source fitting. 

For most of the models we tested using this global fitting method, we treated the 
continuum model as having two components and all of these components are modified by 
Galactic absorption. Since most of the PG QSO spectra do not show strong iron lines (see 
Figure [!])□, we only applied the global model to the continuum. The first of these components 
is AGN emission which includes the standard power law for an AGN spectrum. The second 
is the soft-excess component. This can be a blackbody, a Comptonization model, or a 
reflection model. For all the objects, we linked the photon index, soft-excess temperature, 
and the normalizations for each component to be the same for all objects (i.e. the median 
model for the class). However, we allowed model parameters that describe characteristics 
that may vary in individual sources, such as disk inclination, ionization parameter, and 
intrinsic source absorption, to be free. Since each source also has different brightnesses, a 
multiplicative factor was included for each of the components to adjust for the differences 
in intensity amongst the sources as well as the relative contributions between the different 
model components. Simply, the global model can be described in equation form: Modelj = 
AbsGai,j x Absi ntj j x [Xj x AGN + Yj x SE], where Xj is the multiplicative factor for the AGN 
component that varies depending on the source i and Yj is a similar factor for the soft excess 
component. 



4.2. The Favored Model 

Table summarizes the 13 spectral models we examined as possible universal models 
for the PG QSOs. We first tested a simple power law model for the global QSO fit, model 
A of Table [9j It is a very poor fit to the data, giving a reduced \ 2 ~8.0. We then tested the 
absorpti on-based m o dels f or the origin of the soft excess, models B-E. A scattering model 



(B), see iTeng et al.l (120091 ) for a description, provided a better fit (xt ~3.0) but the model 
offered a poor description of the data at higher energies. A single partial covering absorption 
model (C) also gave a very poor fit (Xv ~3.4). An additional partial covering absorber 
(D) provided a much better fit (xt ~l-8); however, the model severely underestimated the 
spectral flux above ~5 keV. In this model, the nominal photon index of the spectrum was 
~2.8, much steeper than the generally observed range for AGNs (1.6-2.2). By adding a 



5 The lack of strong iron lines in PG QSOs may be due to the X-ray Bal dwin Effect where the inten sity 
of the line is inversely correlated with the total luminosity of the quasars (jlwasawa fc Taniguchilll993l ) or 
gravitational smearing (C06). 
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MEKAL component for the soft excess to the double-partial-covering model (E), we derived 
an even better model to the data (xl ~l-5), but the model still underestimated the flux 
above ~5 keV as shown in Figure [7J 



We then considered the Comptonization model of lPorquet et all (model F; 120041 ) for the 
soft excess. This model gives xl ~2.1. While this is a much better fit to the data than the 
simple power law models, large residuals remain, providing a poor fit both at low (<0.5 keV) 
and at high (>4 keV) energies. 

Finally, we tested the reflection-based models for the soft excess, models G-M. The 
reflection-based models are XSPEC models PEXRAV (a neutral reflector), PEXRIV (an 
ionized reflector), and a blurred reflection model (the C06 model). REFLION component in 
the C06 model is more complex than PEXRIV. PEXRIV considers only bound-free transi- 
tions in the reflected spectrum; REFLION also includes the ionization states and transitions 
for O and Fe ions. We first modeled each reflection-based model (models G-I). These at- 
tempts produced very poor results (xl ^3.0) and suggest that reflection cannot explain the 
soft excess. We then added a redshifted blackbody component to each of the reflection mod- 
els for the soft excess. The reduced x 2 value for the model J is ~1.3, a great improvement 
over the previous models. This model describes well the spectra above ~1 keV, but is a poor 
model below this energy (see Figure [7J. The PEXRIV plus blackbody model (model K) 
appeared to be a good fit to the data (xl ~l-2; Figure UJ), except for an absorption feature 
around 0.7 keV (the atomic transition of O VII or O VIII). Lastly, we evaluated the C06 
model plus a blackbody (model L). Surprisingly, this more complex model for an ionized 
reflector gives a worse fit than the simpler model K (xl ~1.4; Figure UJ). In particular, this 
model is a very poor description of the spectra above ~4 keV. 

Based on the statistics and the residuals of the models, the multi-source fitting method 
indicates that a reflection-based model is favored as the universal model. However, there 
is still a requirement for a blackbody component in order to model the soft excess. One 
reason why the PEXRIV component offers a better fit than the REFLION component may 
be because many of the spectra do not require ionized reflection. Only 10/22 objects (45%) 
require the ionization parameter to be above 30 erg cm s" 1 , the minimum value for this 
parameter in REFLION. The PEXRIV component allows a full range of ionization param- 
eters, starting with a minimum value of erg cm s' 1 . It should be noted that model K 
with the PEXRIV component requires a steeper photon index (2.37) than model L with the 
REFLION component (2.02). While a photon index of ~2.4 is steeper than the generally 
accepted value, it is the median value seen in PG QSOs (see Figure E]). When a redshifted 
absorption edge centered at around ~0.68 keV was added to model K, the statistics greatly 
improved (xl ~1.1; Figured]). For the three objects with EPIC-MOS data that were not 
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included in the multi-source fitting, we applied the best-fit global model with the AGN and 
soft excess components fixed to individually model the MOS spectra. The best global fit of 
model M implies that F = 2.37 ± 0.02 and kT = 0.127 ± 0.001 keV with an absorption edge 
at 0.68 ±0.01 keV and optical depth of 0.31±g;g|. The nominal energy of the absorption edge 
is consistent with the atomic transitions of O VII and O VIII. The reflected portion of the 
AGN contribution is ~6-37%, with a median value of 18%, of the power law luminosity for 
all 25 objects^. Some of the other parameters derived from model M for the PG QSOs are 
listed in Table [TDJ 



4-2.1. Disk Inclination 

C06 suggested that the XMM-Newton data are sufficient to robustly measure the in- 
clination angle of the accretion disk. Neither of the distributions of inclination angles from 
the global PEXRIV nor the constrained C06 models from the same 22 EPIC-pn spectra 
matches a random distribution. The inclination measurements are highly dependent on the 
model, which may indicate that the data cannot adequately assess the inclination values. 
This may also be a result of the small number of quasars in the sample. In any case, we 
caution against relying too heavily on the inclination measurements derived from the mod- 
els whether fitted to individual or multiple sources. Since the subsequent analysis is based 
mainly on good X-ray flux measurements, we choose the model M as the best-fit model 
despite its shortcomings. 



4-2.2. Origin of the Soft Excess 

The global modeling of the quasar spectra suggests that the soft excess is not well 
described by an absorption model (§ I4.2p . Figure M clearly shows that the 0.5-10 keV flux 
from the soft excess component is correlated with that from the AGN component based on 
model M. Linear regression analysis suggests that the two components are linearly correlated: 



6 A reflection percentage of 37 is rather high. The two sources that have 37% reflected fraction are 
PG 1001+054 and PG 1004+130. The former object is the quasar with the fewest detected counts by EPIC- 
pn. Its lack of counts, particularly in the hard band where the reflected component dominates, may have 
resulted in a poor fit. The PG 1004+130 spectrum used in the universal fit is the average of two epochs. As 
can be seen in Figure [T] and unlike other objects presented in this paper with multiple observations, there 
is noticeable change in the spectral shape, particularly at below 2 k e V. Th e average spectrum may not be 
representative of the source spectrum in general. Also, I Miller et al.l (|2006r ) studied this object extensively. 
There is no obvious soft excess in this object. 
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logL bbody = (1.04±0.28)logL p i-(2.69±12.28). The correlation is very strong with R 2 = 0.65 
and a significance of > 99.99%. The outlier, PG 0050+124 (I Zw 1), is removed from the 
regression analysis; recall that this is the only source that retains an absorption feature in 
the residuals after the inclusion of the 0.7 keV absorption edge (see Figure E]). To within 
the errors, the expression derived from the PG QSO data is consistent with tha t found in 
a sample of Swift /BAT AGNs with soft excess (R 2 = 0.48; Iwinter et aD l2009h . but with 
a greater correlation coefficient. PG 1501+106 is the only overlap between the BAT and 
our PG QSO samples. The soft excess in AGNs and PG QSOs seems to arise from the 
same process. The linear relationship between the blackbody and p ower law luminositie s 
precludes absorption as the origin of the soft excess, in agreement with lWinter et al.l ( 120091 ). 
A starburst origin of the soft excess in these objects is also ruled out. The soft excess 
luminosities are much higher than the expected starburst X-ray luminosities based on FIR 
measurements of star-forming galaxies ( IPersic et al.l I20Q4J ) . 



Since the global model fits the AGN and the soft excess components independent of 
each other, it does not require the AGN components to be correlated with the soft excess. 
Therefore, the linear relationship between the input power law and the soft excess luminosi- 
ties is significant and implies a link between the source of the power-law emission and the 
soft excess. We also compared the soft excess luminosity with black hole mass and ionization 
state of the reflector. Neither of these quantities is apparently correlated with the soft excess. 



Many authors (e.g.. iGierlinski &: Dondl2004j ; IWinter et al.ll2009l ) have argued that if the 



soft excess were due to thermal emission arising from the accretion disk, then the blackbody 
temperature should correlate with the mass of the black hole or the Eddington ratio. Indeed, 
in this scenario, the thermal temperature should scale according to T oc M~ 1 ' i (L/LEdd) 1 ' 4: , 
assuming that all of the gravitational energy gained from accretion is re-radiated by the disk. 
Contrary to this, previous studies and our analysis in § 13.11 have found that the thermal 
temperature of the soft excess is consistently ~0.1 keV. The constancy of the observed tem- 
perature is unexpecte d since the masses of both Seyferts (jWinter et al.l 120091 ) and PG QSOs 



( IVeilleux et al.ll2009bl ) span two orders of magnitude (~ 10 7 — 10 9 M Q ). Using the absorption 



corrected 2-10 keV to Eddington luminosity ratio as a proxy for the Eddington ratio, we 
find that this quantity also spans two orders of magnitude for both PG QSOs and Seyferts 
(Figure and lWinter et al.l 120091 ) . Figure [TO] is a histogram of the predicted PG QSO disk 



temper atures calculated using the photometric black hole mass estimates of IVeilleux et al. 



( 12009bl ) and the 2-10 keV luminosities. These are only estimates because the hard X-ray to 
bolometric luminosity correction for AGNs is un certain (e.g., extre me ultraviolet measure- 
ments are not available for most of these objects). iLusso et al.l ( 120101 ) found in 150 COSMOS 
AGNs that the 2-10 keV to bolometric luminosity correction can range from about 10 to 
about 200 with some dependence on the Eddington ratio. For our sources, given the weak 
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dependence of the temperature on the Eddington ratio, the predicted disk temperatures may 
be underestimated by a factor of a few at most. As the figure shows, the disk temperature 
is surprisingly constant for all objects, with a mean of ~1 eV and a standard deviation of 
~0.4 eV. Only a fraction of all QSOs could be plotted in thi s figure due to the lac k of black 



hole mass information for many objects in our sample (see IVeilleux et al.ll2009bf ). so small 
number statistics may be an issue here. Nevertheless, the predicted temperature of the disk 
seems too low to explain the measured soft excess temperature of ~0.1 keV. 

The mechanism producing the soft excess must be mass independent. For example, 
suppose the corona is geometrically thick and thus at a temperature comparable to the 
local virial temperature in all sources. Then T oc y/M/R, which at a given number of 
gravitational radii or a given orbital speed (e.g., the speed characteristic of the broad line 
region) is the same for any mass, hence Compton up scattering could produce a characteristic 
spectrum that is independent of mass. Therefore, coronal effects may contribute to the soft 
excess. However, the main argument against this scenario is that the virial temperatures 
corresponding to broad line region velocities (a few xlO 3 km s -1 ) are much too high (a few 
keV) to produce the soft excess. 

The multi-source fitting has demonstrated that all the models discussed have difficulties 
explaining all the properties of the soft excess emission. 



5. The Multi-wavelength Properties of the QUEST ULIRGs and PG QSOs 

5.1. FIR Classes of PG QSOs 



Netzer et al.l (120071 ) separated PG QSOs into two different FIR classes, strong and weak 
FIR emitters, based on the 60-to-15 /zm luminosity ratio. The authors found that both types 
of QSOs have similar underlying AGN spectral energy distributions (SEDs) in the infrared 
(> 1 jum). We find that this statement seems also valid in the X-rays: The appears to be 
no correlation with the FIR strength with the AGN contribution to the infrared luminosity 
(Figure H]). The C06 models suggest that the strong FIR-emitters have softer spectra (due 
to the possible contribution from a starburst to the 0.5-2 keV band) than the weak FIR- 
emitters, but these models are not favored for the reasons discussed in § 14.21 
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5.2. AGN Contribution to the Bolometric Luminosity 

V09a presented six independent mid-infrared methods for determining the AGN con- 
tribution to the overall bolometric luminosity of ULIRGs and QSOs as part of the Spitzer- 
QUEST study. The details of the methods are presented in the Appendix of V09a and 
not repeated here. Since the hard X-ray emission is dominated by the AGN (assuming no 
Compton-thick sources), we use the ratio of the absorption-corrected 2-10 keV luminosity to 
the bolometric luminosity as a proxy for the AGN contribution to the bolometric luminosity 
and compare the results with those derived by V09a. The results on the QSOs and ULIRGs 
are presented in Figure [TT1 

Figure [11] shows a positive correlation between the Spitzer-derived AGN contribution 
and the absorption-corrected 2-10 keV to bolometric luminosity ratio, although considerable 
scatter is evident in the relation. Note also that the range of hard X-ray to bolometric lumi- 
nosity ratios is considerably broader than the range of Spitzer-derived AGN contributions. 
The Spitzer results are likely more uncertain in objects with intermediate and high r e // (as 
discussed in V09a) , but Figure [11] shows no link between r e f / and the discrepancy between 
the Spitzer and X-ray results. We discuss two possible causes for this discrepancy, aside from 
the obvious and likely possibility of real intrinsic var iations in the hard X-ray to bolom etric 



luminosity ratio among pure AGN (see Figure 11 of iLusso et al.ll2010t I Just et aJj 120071 ). 



5.2.1. Zero-point Calibrations 



One potential source of uncertainty is the definitions of the zero-points for the mid- 
infrared methods in V09a. The pure- AGN zero points were derived from the FIR-undetected 
PG QSOs. Following iNetzer et all (120071 ). the lack o f FIR detection was assumed to indicate 
a lack of starburst in these objects. As discussed in INetzer et al.1 (120071 ). this assumption is 
likely accurate to ±10-30% and thus cannot account for the substantially larger discrepancies 
between the X-ray and mid-infrared methods. The pure-starburst zero points were derived 
from HII ULIRGs, which are known to be different in terms of gas density and radiation 
fields from less luminous, optically-selected starbursts. For ins tance, the iW/iR m icron 
flux ratios and the 7.7 micron PAH equivalent widths derived by iBrandl et al.1 ( 120061 ) for a 
sample of optically-s elected starbu r st nuc lei are lower than those found in the HII ULIRGs. 
If we were to use the IBrandl et al.1 ( 120061 ) values as the pure-starburst zero points, then the 
estimated AGN contributions from the mid-infrared data would be systematically reduced, 
but this still could not explain the very broad range in hard X-ray to bolometric luminosity 
ratio observed among ULIRGs and QSOs. 
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5.2.2. A Matter of Obscuration 

A more likely explanation for the discrepancy between the two methods is related to 
the fact that the hard X-ray and mid-infrared observations probe very different regions of 
the AGN. The hard X-rays represent direct or reflected emission from the accretion disk of 
the black hole or from material near it. The X-ray source occupies a small volume, on the 
scale of less than a parsec. In contrast, the mid-infrared diagnostic method of measuring 
AGN emission relies on the detection of high-ionization fine-structure line or PAH emission 
which is produced from a larger volume (tens to hundreds of parsecs). Thus, the mid- 
infrared emitting region is less likely to be affected by obscuration than the central X-ray 
emitting region. The broad range in hard X-ray to bolometric luminosity ratio may be due 
to unsuspected obscuration that is not correlated with the mid-infrared extinction measured 
on a larger scale (Figures H] and [TT]) . Both Chandra and XMM-Newton operate at 0.5- 
10 keV, below the peak of the Compton reflection hump, the dete ction of whi c h can better 



constrain the absorbing column. Mrk 273 is a good case in point. iTeng et al.l (120090 had a 
marginal detection of Mrk 273 above 10 keV with Suzaku. The simultaneous modeling of the 
Suzaku, Chandra, and XMM-Newton spectra found that the source is highly obscured. The 
absorption corrected 2-10 keV luminosity with the Suzaku data is ~ 3.2 times that derived 
from Chandra or XMM-Newton data alone. Mrk 273 is also one of the sources that exhibits 
moderate extinction from the mid-infrared data (T e // ~ 6.4). This is a change of about 
half an order of magnitude or ~ 30% in AGN percentage (Figure H]). Mrk 231 is another , 



more extreme, example: The detection of hard X-rays above 10 keV by [Braito et al.l (120041 ) 
boosted the absorption-corrected hard X-ray luminosity by more than an order of magnitude, 
bringing it to a value consistent with that of powerful QSOs. In summary, the absorption- 
corrected hard X-ray luminosities of ULIRGs and, to a lesser extent, QSOs are uncertain, 
so it is not inconceivable that the very broad range of hard X-ray to bolometric luminosity 
ratios of ULIRGs and QSOs is due at least in part to unsuspected obscuration of the hard 
X-ray source. 



5.3. Trends with Merger Phase 



While the hard X-ray to bolometric luminosity ratio may not be an accurate absolute 
measure of the AGN contribution to the bolometric luminosity of ULIRGs and QSOs, trends 
of this ratio within these classes of objects may provide useful relative information. 



In the evolutionary scenario of lSanders et al.l (119881 ) . starburst-dominated "cool" U/LIRGs 
evolve into AGN-dominated "warm" ULIRGs and then eventually optically selected quasars. 
In this scenario, the AGNs turn on only near the end of the merging process. The recent 
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Spitzer results of V09a have brought some support to this scenario. The X-ray data also 
seem to be largely consistent with this picture. The top panel of Figure [12] supports previ- 
ous arguments that "warm" ("cool") objects like Seyfert 1 ULIRGs and PG QSOs (HII and 
LINER ULIRGs) are AGN-dominated (starburst-dominated). ULIRGs with intermediate 
infrared colors also have intermediate hard X-ray to bolometric luminosity ratios. This sup- 
ports the idea that the cool, starburst-dominated, and obscured objects evolve to become 
warm, AGN-dominated, "naked" quasars. 

The middle panels of Figure d2] demonstrate that the more X-ray luminous AGN are 
more likely to be sources optically classified as either a QSO or a Seyfert galaxy. The HII 
ULIRGs have the lowest hard X-ray to bolometric luminosity ratios on average, as expected. 
The LINER ULIRGs show the broadest range of 2-10 keV to bolometric luminosity ratios, 
suggesting tha t the energy source of infrared-selected LINERs is a wide mixture of AGN 



and starburst (ISturm et al.ll2006l : IVeilleux et al.ll2009ai ). The very low L2— iokev/ ^boi in some 



ULIRGs may be a sign of unsuspected large X-ray absorptions. 



Veilleux et al.l ( 120021 ) defined interaction classes (IC) for merging systems where class I 
are systems on their first approach. The galaxies are unperturbed, with no evidence of tidal 
features. The IC II represents sources that are in the first-contact phase of the interaction 
where strong bars or tidal tails have yet to be formed. The pre-merger stage, or IC III, 
consists of systems that show identifiable binary nuclei and strong tidal features such as tails 
and bridges. This class is further divided into two sub-categories: IC Ilia are sources that are 
wide binaries where the nuclei are > 10 kpc (projected) apart and IC Illb are sources that are 
close binaries where the nuclei are < 10 kpc (projected) apart. IC IV represents the merger 
stage of the interaction. These sources have prominent tidal tails. The sub-category of IVa 
(diffuse mergers) includes sources that have diffuse central regions whereas IVb (compact 
mergers) includes sources that are dominated by a single nucleus. The division between IVa 
and IVb represents the difference between binary and single nucleus objects. Finally, IC 
V is synonymous with old mergers. These sources do not have strong tidal features, but 
their central morphologies are disturbed, similar to those of objects in IC IV. Therefore, the 
sequence of IC I-V represents the complete merging sequence of galaxy systems. It is of note 
that almost all 1-Jy ULIRGs are classified as IC III-V (moderate age mergers) and the PG 
QSOs as IC IVb-V (old mergers). 

The bottom panels of Figure [12] suggest that AGN activity is often most dominant in 
coalesced remnants corresponding to the latest stage of the merger - IVb and V. The steep 
increase of AGN dominance between classes Illa/b (early stages of the merger) and IV/V 
suggests that black hole growth picks up in the post-merger phase of the interaction, as 
found in the Spitzer data. The scatter among each interaction class indicates that stochastic 
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events may trigger nuclear activity at any time along the merger sequence. However, the 
likelihood of finding AGN-dominated U/LIRGs increases along the merger sequence and in 
objects with warmer infrared colors (/2s//6o)- 

It is inter esting to compare the se results with the predictions of recent numerical sim- 



ulations (e.g., iHopkins et all 120081 . and references therein). These models suggest that 



the starburst domin ates the total luminosity prior to and during the merger (phase D of 
Hopkins et all 120081 ) . After coalescence, the central black hole in these simulations grows 
rapidly before the "blowout" phase (phase E) where an AGN-driven wind is purported to 
expel the remaining dust and gas, removing material for both accretion and star formation. 
The result is a luminous, blue quasar with little star formation (phase F). The accretion rate 
of the active nucleus is predicted to peak between phases D and E. Then the luminosity of 
the quasar fades during the post-blowout quasar stage. For comparison, Figure [13] shows 
the evolution of the AGN luminosity along the final stages of the merger sequence derived 
from the mid-infrared and X-ray data. The data indicate that star formation peaks prior to 
the coalescence of the nuclei, while the accretion rate onto the black hole increases rapidly 
to peak at coalescence. There is no significant change in the quasar luminosity after coales- 
cence. The fading of the quasar must happen after the epoch covered by the QUEST sample 
of QSOs. The purport ed galactic-scale win ds in these simulation s are observed in several 



QUEST ULIRGs (e.g., iRupke et al.ll2005al Jbl; IVeilleux et al.ll2005l and references therein). 



6. Summary 

We have performed an uniform analysis of X-ray data on 40 ULIRGs and 26 PG QSOs 
from the QUEST sample. New and archival observations obtained using Chandra and XMM- 
Newton were used for this study. The X-ray results were compared with those recently 
derived from Spitzer IRS spectra. The major conclusions are as follows: 

1. By fitting the PG QSO spectra simultaneously, we favor a reflection-based model with T 
= 2.37 ± 0.02 and a black-body component with kT = 0.127± 0.001 keV to account for 
the ubiquitous soft excess. The best-fit universal model indicates that the soft excess 
luminosity /flux is linearly related to the 0.5-10 keV absorption-corrected power-law 
luminosity /flux. This implies that the source of the soft excess is directly related to 
accretion onto the central black hole, rather than an external factor such as intervening 
absorption. An absorption edge at 0.68 ± 0.01 keV with optical depth of 0.311q q3 is 
required for the universal spectral model to fit the QSO spectra as a class. This edge 
is consistent with the atomic transitions of O VII and O VIII. 
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2. There does not appear to be any correlation between the FIR emission strength of the 
QSOs and th e 2-10 keV to bolometric luminosity ratio. This extends to the X-rays the 



conclusion of iNetzer et al.l ( 120071 ) that there is no obvious difference in the underlying 



infrared AGN SED of strong and weak FIR emitting PG QSOs. 

3. Using the absorption-corrected hard X-ray luminosity as a proxy for the AGN con- 
tribution to the bolometric luminosity in these systems, we find results that generally 
agree qualitatively with those from V09a. The ratio of absorption-corrected hard X-ray 
to bolometric luminosity is not an accurate absolute measure of the AGN contribution 
to the bolometric luminosity of ULIRGs and QSOs, but trends of this ratio within 
these classes of objects provide useful relative information. The likelihood of powerful 
nuclear activity increases along the merger sequence and in objects with warmer dust 
temperatures and AGN optical signatures. The scatter in these trends is likely due to 
stochastic accretion events, unsuspected X-ray absorption, or variations in the intrinsic 
hard X-ray to bolometric luminosity ratio of pure AGNs. 

4. The bolometric luminosity of the AGN in U/LIRGs and PG QSOs evolve with merger 
stage. The starburst seems to dominate the total power prior to the merger. Then 
the accretion rate onto the black hole increases rapidly during coalescence at which 
point the AGN dominates the bolometric luminosity. The predictions from numerical 
simulations are largely consistent with these results. 
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APPENDIX A — The Reliability of the Hardness Ratio Method 

As discussed in § 13.2} U/LIRGs are notoriously difficult to observe in the X-ray. Most of 
these objects are faint, either due to the lack of an AGN or the presence of heavy obscuration. 
Traditional spectral fitting cannot be used to model the complex s pectra of th e se sou rces 



when counts are limited. The hardness ratio method developed in iTeng et al.l (120051 ) has 
proven to be effective in finding obscured AGNs in at least one case. IRAS F04103-2838 
was found to contain both a starburst and an AGN with only 30 detected counts in a 10 ksec 
Chandra observation. A deeper (~20 ksec) XMM-Newton follow-up revealed an Fe K a line 



at rest-frame energy of ~6.5 keV, consistent with cold neutral iron f lTeng et al.ll2008f ). The 
hard X-ray emission is dominated by a nearly Compton-thick AGN with intrinsic 0.2-10 keV 
luminosity ~10 43-44 ergs s _1 . 

Since the detected counts of the majority of objects in the U/LIRG surveys are low, the 
errors associated with the HR method are inherently large. To further test the reliability of 
the HR method in recovering the input spectrum, we have performed a set of simulations. In 
these simulations, we set the input model of an unabsorbed AGN at z ~0.1 as a redshifted 
power law with T at 1.8 (the canonical value for AGNs) absorbed by a Galactic column of 
2 x 10 20 cm -2 . The normalization of the input power law model was such that the model 
0.5-10 keV flux is ~ 5 x 10~ 14 ergs s _1 cm -2 , the mean value observed in faint U/LIRGs. 

We first tested the dependence of the HR method on the number of detected counts by 
varying the exposure times in the simulations. For each exposure time tested, the average of 
1000 simulations showed that the nominal T and the 0.5-2 keV and 2-10 keV flux values are 
remarkably stable even when the "detected" counts were reduced to as low as 30 (a 5 ksec 
exposure). On average, the output T determined from the HR method remained the same 
as the input, but the error bars increased as the number of detected counts decreased. The 
nominal 0.5-2 and 2-10 keV fluxes are within 1% of those from the input model. The top 
portion of Table [TT1 summarizes these results. 

In the HR method, we assumed that the only absorption is from the Galaxy, but this is 
not the case for the U/LIRGs we observe. The internal absorption in these objects are often 
high. Therefore, we also tested the dependence of the HR method on the intrinsic absorption 
of the source. To this end, we added intrinsic source absorption to the input spectrum. For 
15 ksec exposures, we varied the intrinsic source absorption between io 20-22 cm -2 by steps 
of 5 x 10 20 cm -2 to see how the measured spectral parameters using the HR method are 
affected. Figure dH shows the results of these simulations. As the internal column increases, 
the photon index becomes flatter and then inverted (r < 1). T begins to become inverted 
when the column is > 5 x 10 21 cm -2 as seen in the bottom portion of Table [TTJ Since the 
0.5-2 keV flux is more readily affected by absorption, it deviates from the input spectrum 
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by more than 10% when the absorption is at only 1 x 10 21 cm~ 2 . On the other hand, 
the 2-10 keV flux is more stable, but begins to deviate from the input spectrum when the 
absorption is high enough that T becomes inverted. Of course, if the source has a softer 
spectrum (e.g., containing a commonly seen soft-excess component), then it would require a 
higher column for the spectrum to become inverted. 

One of the important usages of the HR method is to estimate the photon index of the 
X-ray spectrum which in turn gives estimates of the 0.5-2 and 2-10 keV fluxes of our targets. 
The spectral index is the parameter that is the basis for the flux estimates. As shown by the 
simulations, the accuracy of T based on the HR method depends on the intrinsic absorption 
of the source. Figure [15] compares the T derived from the HR method and the traditional 
method for the moderately bright U/LIRGs in § 13.21 and in the archived Chandra sample 
discussed in Appendix B. The figure shows that our measurements from the two methods 
are consistent with each other unless the intrinsic column densities are > 10 22 cm -2 , in 
agreement with the simulations. The next consideration is to see how well the estimated 
fluxes from the HR method match those of the traditional method. We have plotted in 
Figure [TBI the 0.5-2 keV and 2-10 keV fluxes derived from both the HR and the traditional 
fitting method with more complex spectral models for the moderately bright U/LIRGs in 
§ 13.21 As the figure shows, the majority of the HR flux values from a simple unabsorbed 
power law model are within 50% of the spectral fitting values of more complex models. The 
HR method is more likely to overestimate the fluxes in both the soft and hard bands. The 
medians values for FHR/Ff its are approximately 1.2 and 1.3 for the soft and hard bands, 
respectively. 



As lTeng et al.l ( 120091 ) and other authors have demonstrated, the X-ray spectra of U/LIRGs 
are often more complex than a simple power law model can characterize. We further tested 
the HR method by assuming that the input model contains a reflection component or that 
the intrinsic absorption is due to a partial covering absorber. The addition of a reflection or 
a partial covering component flattens the AGN power law spectrum, mimicking the effects 
of high intrinsic absorption. For an intrinsic absorption of 1 x 10 22 cm -2 , the HR method 
recovers the input flux very well for both the reflection and partial covering models^. In 
fact, for both models, the recovered 0.5-2 keV fluxes are within 2% of the input model flux. 



7 The reflection model assumes a power law plus a PEXRAV component. The normalization of the 
reflected component is set to 2% of the intrinsic p ower law component as measured in the complex Suzaku 
spectrum of IRAS F05189-2524 (|Teng et al.ll2009l) . The photon index of the power law spectrum is fixed at 
1.8 and the recovered T is 0.701"' Jg. For the partial co vering: model, the covering factor is assumed to be 
90%, similar to that found in Mrk 273 (jTeng et al.ll2009[ ). In this case, the input T is again fixed at 1.8 and 



the HR estimate of V is O^l^g. 
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On the other hand, the HR estimates of the 2-10 keV fluxes are above the input values by 
56.4 and 35.3% for the reflection and partial covering models, respectively. These errors are 
comparable to or better than the estimates obtained for the simple power law model at the 
same column density (Table [TTT) . Therefore, even for objects with complex spectra, the HR 
method is able to provide fair approximations of their spectral properties. 
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Table 1. The Sample 
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Table 1 — Continued 



Galaxy 
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References. — 1 = Chandra AO 10 object; 2 = Teng et al. 2005; 3 = Piconcelli 
et al. 2005; 4 = Crummy et al. 2006; 5 = Teng et al. 2008; 6 = Ptak et al. 2003; 7 
= Imanishi 2004; 8 = Teng et al. 2009; 9 = Franceschini et al. 2003; 10 = Sanders 
et al. (2003); 11 = XMM-Newton AO 7 object 

Note. — Col.(l): Galaxy name. Coordinate- based names beginning with "F" are 
sources in the IRAS Faint Source Catalog. Col. (2): Redshift. Col. (3): Bolometric 
luminosity. For ULIRGs, we assume L(bol) = 1.15L(IR ). For PG QSOs, we assume 
L(bol) = 7L(5100 A) + L(IR) fsee lNetzer et alj|2007t for details on L(5100 A)). 
Col. (4): Optical spectral type, from Veilleux et al. (1995, 1999a) and Rupke et al. 
(2005a). Col.(5): Interaction class, from Veilleux et al. 2002, Veilleux et al. 2006, 
or Veilleux et al. 2009 (in order of preference). Col. (6): Nuclear separation, in kpc. 
Col. (7): Galactic column density from Dickey and Lockman (1990) in units of 10 20 
cm -2 . Col. (8): References for X-ray data. 
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Table 2. Journal of Observations 
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ZUU1 IVldilOIl _L / 


YA/TA/T 


00^^6401 01 
UUOOU"4:U1U1 


Hcckm9.n 


6 3 




9001 Ortnhpr 30 

ZUU.L V/LliUULl UU 


exo 


6001 6S9034 


r— 1 /"^ c lr TYi '"i 
1 L Lis. 1 1 1 CL 1 1 


19.3 




9009 Tannarv 30 


LAU 


6001 683439 


HcckrQ9.ii 


14 7 


F07 I ^Q84-6 C >08 


9000 A/Tarrh 91 
zuuu ivxcu.1^11 z J- 


exo 


7001 91 K1 6 

/ UU J- Zr J.O ±U 


Green 


1.5 




9001 Ortnhpr 9^ 


YA/TA/T 


00Q4400301 


oancicrs 


19 7 


F08^79-U3Q1 ^ 


9004 Anvil 1 3 


XMM 


09006301 01 

UZUUUOU ±U -L 


T t > i aTiiohi 
.LlllellllSIll 


13.3 




9006 Tannarv 96 

ZiUUU -J till Llcti y ZiU 


exo 


701 9K06S69 

/ U-LZOUUOUZ 


I < , 1 1 [ 1 fifino 


14.3 


T?nonQO i nc;nQ 
r UyUoy+UoUo 


zUUo ucccmbcr ol 




< UlooyiUo4D 


Vcillcux 


1 fi o 


UGC 5101 


2001 May 28 


exo 


6001672033 


Heckman 


48.7 




2001 November 12 


XMM 


0085640201 


Hcckman 


24.0 


F09539+0857 


2004 January 6 


exo 


7008894806 


Vignali 


5.2 


F10190+1322 


2003 January 31 


exo 


7007864107 


Wilson 


9.4 




2003 May 5 


XMM 


0146990101 


Risaliti 


16.3 


F10565+2448 


2003 October 23 


exo 


6003263952 


Martin 


29.2 




2003 June 17 


XMM 


0150320201 


Martin 


22.9 


F11095-0238 


2009 April 9 


exo 


70189010347 


Vcillcux 


15.5 


F11119+3257 


2002 June 30 


exo 


7005763137 


Anabuki 


15.2 


F11223-1244 


2006 June 8 


XMM 


0405950101 


Terashima 


0.0 


F12072-0444 


2003 February 1 


exo 


7007884109 


Wilson 


9.2 


F12112+0305 


2001 December 30 


XMM 


0081340801 


Franccschini 


16.0 




2003 April 15 


exo 


7007894110 


Wilson 


10.0 


3C 273 a 


2000 June 14 


XMM 


0126700301 


Janscn 


44.7 




2000 June 17 


XMM 


0126700801 


Janscn 


36.0 




2001 June 13 


XMM 


0136550101 


Janscn 


62.0 




2003 July 7 


XMM 


0159960101 


Pollock 


40.6 




2004 June 30 


XMM 


0136550801 


Janscn 


11.7 




2007 January 12 


XMM 


0414190101 


Parmar 


47.9 


Mrk 231 


2000 October 19 


exo 


7003891031 


Garmire 


36.0 




2001 June 7 


XMM 


0081340201 


Franccschini 


17.2 
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Table 2 — Continued 



Galaxy 


Date 


Telescope 


Obs ID 


PI 


Exposure 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Mrk 231 (cont.) 


2003 February 3 


exo 


7007074028 


Gallagher 


41.0 




2003 February 11 


exo 


7007084029 


Gallagher 


49.5 




2003 February 20 


exo 


7007094030 


Gallagher 


32.2 


F13218+0552 


2004 July 11 


XMM 


0200660301 


Matt 


5.1 


Mrk 273 


2000 April 19 


CXO 


700114809 


Xia 


41.0 




2002 May 7 


XMM 


0101640401 


Aschenbach 


18.0 


F13451+1232 


2000 February 24 


CXO 


700141836 


Ward 


19.7 


F14348-1447 


2002 July 29 


XMM 


0081341401 


Franccschini 


6.8 




2006 March 12 


CXO 


7012796861 


Komossa 


14.8 


F15130-1958 


2003 June 2 


CXO 


7007904111 


Wilson 


9.8 


F15250+3608 


2002 February 22 


XMM 


0081341101 


Franccschini 


15.0 




2003 August 27 


CXO 


7007914112 


Wilson 


9.2 


Arp 220 


2000 June 24 


CXO 


700174869 


Clements 


56.1 




2002 August 11 


XMM 


0101640801 


Aschenbach 


10.5 




2003 January 15 


XMM 


0101640901 


Aschenbach 


8.3 




2005 January 14 


XMM 


0205510201 


Sanders 


5.8 




2005 February 20 


XMM 


0205510401 


Sanders 


4.8 




2005 February 27 


XMM 


0205510501 


Sanders 


0.0 


F15462-0450 


2009 April 23 


CXO 


70189110348 


Veilleux 


15.2 


F16090-0139 


2003 February 10 


CXO 


7007924113 


Wilson 


9.8 


NGC 6240 


2000 September 22 


XMM 


0101640101 


Aschenbach 


11.3 




2001 July 29 


CXO 


7002061590 


Predehl 


37.1 




2002 March 12 


XMM 


0101640601 


Aschenbach 


5.8 




2003 March 14 


XMM 


0147420201 


Netzer 


3.4 




2003 March 18 


XMM 


0147420301 


Nctzcr 


0.0 




2003 August 13 


XMM 


0147420401 


Netzer 


7.8 




2003 August 21 


XMM 


0147420501 


Nctzcr 


3.4 




2003 August 29 


XMM 


0147420601 


Netzer 


1.6 


F17208-0014 


2001 October 25 


CXO 


6001692035 


Hcckman 


49.0 




2002 February 19 


XMM 


0081340601 


Franccschini 


9.8 




2003 May 7 


CXO 


7007934114 


Wilson 


8.6 


F21208-0519 


2009 March 31 


CXO 


70189210349 


Veilleux 


15.1 


F21329-2346 


2009 June 21 


CXO 


70189310350 


Veilleux 


15.2 


F22491-1808 


2001 May 24 


XMM 


0081340901 


Franccschini 


17.9 




2007 July 13 


CXO 


7014857821 


Sanders 


14.4 


F23234+0946 


2009 August 15 


CXO 


70189410351 


Veilleux 


15.1 


PG QSOs 


PG0050+124 


2002 June 22 


XMM 


0110890301 


Watson 


18.3 




2005 July 18 


XMM 


0300470101 


Gallo 


57.8 


PG0838+770 


2009 March 2 


XMM 


0550870401 


Veilleux 


12.5 


PG0844+349 


2000 November 5 


XMM 


0103660201 


Aschenbach 


8.5 


PG0953+414 


2001 November 22 


XMM 


0111290201 


Watson 


10.9 


PG1001+054 


2003 May 4 


XMM 


0150610101 


Schartcl 


8.7 
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Table 2 — Continued 



Galaxy 




Date 


Telescope 


Obs ID 


PI 


Exposure 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


PG1004+130 


2003 


May 4 


XMM 


0140550601 


Brandt 


18.0 




2005 


January 5 


exo 


7010295606 


Garmirc 


41.3 


PG1116+215 


2001 


December 2 


XMM 


0111290401 


Watson 


5.5 




2004 


December 17 


XMM 


0201940101 


Stccnbruggc 


5.6 




2004 


December 19 


XMM 


0201940201 


Steenbruggc 


5.0 


PG1126-041 


2004 


December 21 


XMM 


0202060201 


Schartcl 


28.7 


PG1211+143 


2001 


June 15 


XMM 


0112610101 


Turner 


48.9 




2004 


June 21 


XMM 


0208020101 


1! ceves 


34.4 


PG1229+204 


2005 


July 9 


XMM 


0301450201 


Matt 


17.2 


PG1244+026 


2001 


June 17 


XMM 


0051760101 


Fabian 


4.3 


PG1307+085 


2002 


June 13 


XMM 


0110950401 


Watson 


10.2 


PG1309+355 


2002 


June 10 


XMM 


0109080201 


Mason 


23.9 


PG1351+640 


2004 


June 23 


XMM 


0205390301 


O'Brien 


43.3 


PG1411+442 


2002 


July 10 


XMM 


0103660101 


Aschcnbach 


21.9 


PG1426+015 


2000 


July 28 


XMM 


0102040501 


Jansen 


0.5 




900Q 


February 2 




0^08070901 


Veilleux 


A 


PG1440+356 


2001 


December 23 


XMM 


0107660201 


Mushotzky 


15.3 




2003 


January 1 


XMM 


0005010101 


Guainazzi 


17.2 




2003 


January 4 


XMM 


0005010201 


Guainazzi 


10.4 




2003 


January 7 


XMM 


0005010301 


Guainazzi 


18.1 


PG1448+273 


2003 


February 8 


XMM 


0152660101 


Kawaguchi 


18.1 


PG1501+106 


2001 


January 13 


XMM 


0112910201 


Turner 


5.9 




2001 


January 13 


XMM 


0070740101 


Petrucci 


7.6 




2001 


January 14 


XMM 


0070740301 


Petrucci 


9.0 




2005 


January 16 


XMM 


0205340201 


Petrucci 


29.7 




2005 


July 17 


XMM 


0205340401 


Petrucci 


16.3 


PG1613+658 


2001 


April 13 


XMM 


0102040601 


Jansen 


1.8 




2001 


August 29 


XMM 


0102041301 


Jansen 


2.1 


PG1626+554 


2002 


May 5 


XMM 


0109081101 


Mason 


0.4 


PG2130+099 


2003 


May 16 


XMM 


0150470701 


Santos-Lleo 


25.4 


B2 2201+31A 


2008 


December 1 


XMM 


0550871001 


Veilleux 


9.3 


PG2214+139 


2002 


May 12 


XMM 


0103660301 


Aschcnbach 


6.6 



Note. — Col.(l): Galaxy name. Coordinate-based names beginning with "F" are sources in 
the IRAS Faint Source Catalog. Col. (2): Observation start date. Col. (3): Telescope used (CXO 
= Chandra; XMM = XMM-Newton). Col. (4): Observation identification number. Col. (5): 
Principal investigator of the observation. Col. (6): Net exposure time in ks after selection of 
good time interval. 

a 3C 273 is often observed by XMM-Newton for calibration purposes. We selected the six 
observations that have the longest integration times and also maximize the period of elapsed 
time between the first and final observations. 
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Table 3. PG QSOs: PL Model Best-fit Values 



Galaxy 


Obs 


. Date 


kT 




r 


N H 


E line 


"line 


EW 


Xl (c 


.o.f.) 


-Fo.5-2 


^2-10 


-^0.5-2 


^2-10 


(!) 




(2) 


(3) 




(4) 


(5) 


( b ) 


(7) 


(8) 


0) 


(10) 


(n) 


(12) 


(13) 


PG0050+124 


2002 


Jun 22 a 




, ,.+0.02 
z " 31 -0.02 


04+ 01 
u - u ^_ 03 








1.22 


630) 


„ ,o+9.31 
8-38_ 07 


„ ,5 + 7.57 
8.25_ 07 


0.96 


0.74 




2005 


Jul 18 a 




, ,,+0.009 


B6 - 86 Il4.59 








1 09 


812) 


, ,,. + 0.03 
b l °-0 03 


4 93+0.11 
^• yo — 17 


0.93 


0.68 


PG0838+770 


2009 


Mar 2 




1 


^ -0.08 










1.06 


138) 


29+ 03 
u - zy _ 04 


n 7s +0.08 
0- '8_o 08 


0.14 


0.35 


PG0844+349 


2000 


Nov 5 


15+ 005 
u ' iJ -0.06 


2 


,,,, + 0.05 
66 -0.06 


6 13+ 3 - 03 


6.4(f) 


0.0(f) 


o-ii±S:i? 




264) 


c- cc +0.13 
6.55_ 02 i 


5 13+0-31 
0,1 -0.36 


1.42 


0.63 
















6.7(f) 


0.0(f) 


0.18±° ; 2 7 














PG0953+414 


2001 


Nov 22 


16+ 003 
u ' lb -0.07 


2 


44+0-03 
— 0.03 








1.04 


324) 


3 39+ 05 
J " :1M -0.07 


2 93+ 17 


10.89 


6.51 


PG1001+054 


2003 


May 4 


09+ 01 
-0.01 


2 


01 + 0.67 


8.09t 5 3 ;f 7 








0.96 


(19) 


02+ 01 
u - uz _ 01 


19+0- 05 
u — 09 


0.10 


0.12 


PG1004+130 


2003 


May 4 




1 


g 7 +0.20 
D '-0.11 


2 99+ 2 ' 67 








1.04 


111) 


09+ 01 
uuy — 0.01 


31 + 0-06 
uol _0.08 


0.30 


0.57 




2005 


Jan 5 




1 


52 +0.17 
° -0.26 


1 44+0-64 
— 0.69 








1.06 


(95) 


o io+°- 01 

— 0.03 


43+ ' 07 
— 0.17 


0.32 


0.77 


PG1116+215 


2001 


Dec 2 


n ns+ - 041 

u - us _0.004 


2 


53 +0.04 
° — 0.03 


27 21+ 16 01 
— 11 26 








1.00 


255) 


4 34+0-22 
0.06 


„ 27+0.23 
° ' —0.29 


8.45 


4.44 




2004 


Dec 17 


„ „„ + 0.007 
— 006 


2 


49+0.01 
-0.01 


01 «1+5.14 
31-61I 4 13 








1.15 


664) 


3 70+0- 03 
°- ' u — 03 


3 H+0-07 

— 09 


8.35 


4.65 




2004 


Dec 19 


12+ ' 13 
u — 04 


2 


+ 0.04 
o± _0.04 


20 21+ 5 ' 94 
zu - z± — 5 10 








0.98 


238) 


4 00+ -' 09 
'■ uu _ 10 


, 5,+0.'23 
°- 0± — 26 


8.51 


4.63 


PG1126-041 


2004 


Dec 21 


11 + "' 01 
u — 0.01 


1 


„5 + 0.10 
y -0.10 


4 66+ ' 42 








1.13 


309) 


09+ - 01 
u - uy _ 0.01 


1 14+0- 14 
— 0.28 


0.11 


0.13 


PG1211+143 


2001 


Jun 15 a 


19+ 001 
u — 0.02 


2 


oo+O.Ol 
° -0.02 


12.98_ .9 








1.41 


489) 


9 KQ + 0-03 

2-oa_ 04 


2 90+ 08 
z - yu _0.07 


2.13 


0.87 




2004 


Jun 21 a 


21 + 02 


2 


„,+0.02 
b -0.02 


12 40+ 1 ' 64 


6.4(f) 


0.2(f) 


0-18±oT 2 


1.21 


533) 


, ,, + 0.04 
°- AO -0.04 


3 „ 7 +0.09 


1.59 


0.78 


PG1229+204 


2005 


Jul 9 


09+ 01 


2 


,o + 0.03 
dS -0.03 




6.4(f) 
6.7(f) 


0.0(f) 
0.0(f) 


06+ 08 
u - uo _0.06 

11+0- 11 
— 0.09 


1 13 


402) 


„ ,_+0.05 
J± '-0.06 


3 16+0- 16 
J,ib -0.17 


0.61 


0.41 


PG1244+026 


2001 


Jun 17 


62+ 06 
u '°^_0.07 


2 


o +0.07 
72 -0.07 










0.89 


240) 


K r „ + 0.21 
b ' 58 -0.20 


2 55 + O.I3 
^• 00 _0.13 


0.38 


0.14 


PG1307+085 


2002 


Jun 13 


11+ 001 
— 0.01 


1 


oo+O.ll 
oy -0.10 










1.13 


145) 


77+O- 04 
— 04 


1 89+ 23 
±.oa_Q 3 7 


1.05 


1.43 


PG1309+355 


2002 


Jun 10 a 




2 


19 + 0.07 
±y -0.06 


6.02t 3 ;^ 


6.4(f) 


0.0(f) 


12+ 017 


1.03 


264) 


40+ - 02 
u -^ u _ 01 


69+ 07 
u.oy_ 08 


0.74 


0.76 


PG1351+640 


2004 


Jun 23 a 


16+ 007 
u - 1Q — 0.05 


2 


42 +0.04 
* -0.04 


14 - 61 t!:s 2 i 








1.06 


281) 


59+ ' 01 
u - oy _ 0.03 


62 + ' 05 
u - oz _0,06 


0.25 


0.17 


PG1411+442 


2002 


Jul 10 


12 + 003 
ulj! -0.02 


2 


41+0-18 
* -0.15 


26.29t 3 ;™ 


6.4(f) 


0.0(f) 


17+0- 32 
u,i '-0.17 


1.01 


139) 


n ns+O- 007 
0.08_ on 


50+ - 07 
u - ou -0.17 


0.40 


0.26 


PG1426+015 


2000 


Jul 28 




2 


,0 + 0.06 
•"-0.06 










1.01 


148) 


o 14 + 0.35 
° -0.38 


5 9 9 + 0-60 
°' y -0.57 


1.65 


1.09 


PG1435-067 


2009 


Feb 2 




2 


og+O.ll 
Jb -0.10 










1.25 


(60) 


2 60+ 23 
^■ DU -0.25 


2 10+ 037 
^ ,iU -0.35 


1.30 


0.88 


PG1440+356 


2001 


Dec 23 




3 


n9 +0.04 
u ^-0. 04 


8.92+|;j; 








1.17 


164) 


5 , 9 +0.16 
°- 1 ^-0.14 


2 41+0- 28 
-0.31 


2.12 


0.55 




2003 


Jan l a 




2 


Sb -0.01 


14.40+. 2 ' 88 








1.34 


433) 


5 49+0-06 
°'* y -0.05 


„ ,,+0.10 
olA -0.10 


2.37 


0.79 




2003 


Jan 4 




2 


98 + 0.02 
yo -0.02 


„ , r +2.03 
8-25-1.49 








0.96 


344) 


. sfi + 0.07 
4 -86_ .o7 


2 46+°- 16 
z '* b -0.15 


1.99 


0.55 




2003 


Jan 7 




2 


SR + 002 
Sb _0.02 


12 67+ 3 45 








1.45 


340) 


2 90+ - 05 
zyu -0.05 


1 74+O.09 
L - '*-0.10 


1.34 


0.44 


PG1448+273 


2003 


Feb 8 


23+ 01 


2 


„„ + 0.02 
SU -0.01 


16 72+ 6 24 
1D ' '^-4.39 








1.22 


(381) 


3 91+0-05 
° ,y -0.06 


1 97+0- 10 
1 ' y '-0.11 


0.89 


0.31 


PG1501+106 


2001 


Jan 13 a 


°' 08 -0.00 


2 


46+0-02 
* -0.02 


23 12+ 6 ' 40 








1.23 


481) 


14 14+ ' 18 
"■"-0.07 


13 24+ ' 25 


1.12 


0.68 




2001 


Jan 13 a 


09+ 01 

u ' uy -0.01 


2 


_„+0.02 
ou -0.02 


18.66+. 3 '™ 


6.4(f) 


0.0(f) 


04+ 06 
u ' u *-0.04 


1.09 


'566) 


16.08_g'?f 


14 76+ 26 


1.23 


0.70 




2001 


Jan 14 a 


09+ 01 
uuy -0.01 


2 


45+O.02 
^ -0.02 


1 5 SS+ 2 -' 42 

15.88_ 2 07 








1.14 


(624) 


16.02lo.i3 


1 c ,,0 + 0.27 
15.63lo.22 


1.17 


0.71 




2005 


Jan 16 a 


„ „„ + 0.001 
u,uy — 0.001 


2 


2 6 +0.01 
^°-0.01 


13 01+ ' 57 
la ' ul -0.50 


6.7(f) 


0.2(f) 


14+ - 07 
u,± -0.06 


1.56 


821) 


4 93+ ' 04 
*' y,:1 -0.05 


9 78+0- 13 
y ' ,o -0.14 


0.55 


0.47 




2005 


Jul 17 a 


11+ - 01 


2 


,1+0.02 
-0.02 


11 4S+ 1 - 89 


6.7(f) 


0.15(f) 


11+0- 07 
u,ii -0.07 


1.24 


691) 


a sn +0.06 
b '8°_0.09 


11 51+ ' 27 
ll.0-o.29 


0.42 


0.45 


PG1613+658 


2001 


Apr 13 




1 


95+0.10 
y -0.10 


98 94+107.78 








0.99 


135) 


9 5O + 0.38 
2-58_ i9 


4 22+ ' 73 
*'^^-0.82 


1.90 


2.38 




2001 


Aug 29 




2 


,, + 0.08 
1 -0.08 


10 45 + 10 78 








0.98 


246) 


4 00+ 27 


5 29+ - 60 
°-^ y _0.73 


2.74 


2.67 


PG1626+554 


2002 


May 5 




2 


„4+0.15 
u *-0.14 










0.95 


(33) 


9 qq+0.57 
^■ yy _0.48 


3 45 + 0.65 
3- 40 _0.53 


1.49 


1.62 


PG2130+099 


2003 


May 16 a 


n ns+ 004 

U - UB _0.000 


2 


,q + 0.05 
/H -0.05 


- ,,,+0.73 
Oyl -0.62 


6.4(f) 


0.1(f) 


15 + 011 
U10 -0.10 


1.31 


481) 


1 88 + 0° S 
1 '88_o.04 


3 43+0-1S 


0.51 


0.42 


B2 2201+31A 


2008 


Dec 1 


19 + 0.03 
uly -0.05 


2 


,,+0.04 
-0.06 


40 42+ 28 98 
'"■"-16.91 






1.04 


212) 


2 55+O.O7 
z ' oo -0.08 


4 61+ 33 
*' D1 -0.43 


20.18 


19.57 


PG2214+139 


2002 


May 12 


17+ - 02 


1 


sn +0.l6 

8U -0.16 


4 48+0-68 
*'*°-0.68 








1.14 


207) 


n O8+0.03 
"■28_ .04 


3 19+ - 41 
J ' iJ -1.04 


0.29 


0.44 



a Thc spectrum requires additional components for rcdshiftcd absorption edges originating from O VII or O VIII. 



Note. — (f) denotes a fixed value. Col.(l): Galaxy name. Col. (2): Observation start date. Col. (3): Gas temperature from the MEKAL component 
in units of kcV. Col. (4) : Photon index of the direct and indirect power-law components. Col. (5) : Hydrogen column density within the source in units 
of 10 22 cm -2 . Col. (6): Rest-frame energy of the iron emission line(s) in units of keV. Col. (7): Width of the iron linc(s) in units of kcV. Col. (8): 
The equivalent width of the iron line(s) in units of keV. Col. (9): Reduced x 2 °f the best-fit model followed by the number of degrees of freedom 
in (). Col. (10): Observed 0.5-2 kcV flux in units of 10~ 12 erg s -1 cm -2 . Col. (11): Observed 2-10 kcV flux in units of 10~ 12 erg s -1 cm -2 . 
Col. (12): Nominal absorption corrected 0.5—2 keV luminosity in units of 10 44 erg s — 1 . Col. (13): Nominal absorption corrected 2—10 keV luminosity 
in units of 10 44 erg s~ 1 . 



Table 4. PG QSOs: P05 Model a Best-fit Values 
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a Thc best-fit models presented in Table 10 of [Piconcclli ct al. (200^) are used as the basis for the modeling of each spectrum. 
^The spectrum requires additional components for redshifted absorption edges originating from O VII or O VIII. 

c Although Piconcclli ct al. (20(35) found that the double blackbody plus power law model best describes the spectrum, we found that a single blackbody plus power law model is a better fit. 
d Thc absorber here is an ionized absorber with ionization parameter at 237^g^ erg cm 2 . 
c Thc absorber here is a partial- cove ring absorber with the covering fraction at 0-£)6^q'q^. 



* Piconcclli ct al. (2005) suggests the best-fit model is a two-phase warm absorber with a power-law component. The column density parameters are listed in order for the respective phases 
at 1.5 X 10° and 3 X 10 4 K. These two phases have ionization parameters of 99_-^y and 6^^ 2 erg cm 2 , respectively. 



Note. — (f) denotes a fixed value. Col.(l): Galaxy name. Col. (2): Observation start date. Col. (3): Tcmpcraturc(s) of the blackbody, brcmsstrahlung, or Raymond-Smith componcnt(s) 
in units of keV. Col. (4): Photon index of the power-law componcnt(s). Col. (5): Hydrogen column density within the source in units of 10 cm - . Col. (6): Rest-frame energy of the iron 
emission linc(s) in units of keV. Col. (7): Width of the iron line(s) in units of keV. Col. (8): The equivalent width of the iron line(s) in units of keV. Col. (9): Reduced x 2 of the best-fit model 
followed by the number of degrees of freedom in (). Col. (10): Observed 0.5—2 keV flux in units of 10 — 12 erg s — 1 cm -2 . Col. (11): Observed 2—10 keV flux in units of 10 12 erg s 1 cm 2 . 
Col. (12): Nominal absorption corrected 0.5—2 keV luminosity in units of 10 44 erg s —1 . Col. (13): Nominal absorption corrected 2—10 keV luminosity in units of 10 44 erg s 1 . 
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42.8+ 4 ; 2 


0.1+° 


7 


109t 4 6 7 


1.02 


1140) 


2 69+ 1 - 12 
Z - Da -1.65 


4 i9+!- 8 ° 
^■ 13 -4.19 


1 


18 


1.84 




2001 


Aug 29 


1 


89+° 


05 
05 




5.6+- 


1 


8+1 


5 
8 


29.9+™ 


0.1+° 


3 


37lt 77 I 


0.96 


(254) 


4 09+ 26 87 
* ut, -4.09 


4 75 + 22.29 


1 


79 


2.08 


PG1626+554 


2002 


May 5 


1 


74+0 


24 
29 

) 




8 9+ 1 ' 1 


1.6+° 


3 
4 


+34:5 


0.1+ 9 


9 


o 1Q o + 6807 
,3 - Ly,3 -3163 


0.97 


(30) 


oa + 642.02 

2.»a_ 288 


, ,, + 207.31 
^•^a_2.29 


1 


35 


1.07 


PG2130+099 


2003 


May 16 a 




1.00(1 




7 S+0.4 


2-0t° 


1 
1 


31.8+11 


10.0n c 


.2 


B796+|l| 


1.41 


(477) 


1 89+ 1108 
loa -1.31 


3 47 +10.04 
-3.47 





18 


0.33 


B2 2201+31A 


2008 


Dec 1 


1 


71+° 


10 
08 




9 8+°' 2 


1 


2 +o 


1 


71 3+ 7 - 3 


0.3+° 


1 
1 


587 -229 


1.03 


(210) 


2 54+ 2 - 40 
^ 0<l -1.96 


4 62+ 2 ' 60 
*- D - i -3.14 


7 


05 


12.82 


PG2214+139 


2002 


May 12 a ' i ' 


2 


71 +0.04 
' -0.06 


„ „o+0.01 
0.38_ Q 01 


10-0-1.1 


2 


2 +b 

Z -0 


1 

2 


27.4+i'J 


0-1+° 


2 


16lt 2 17 


1.05 


(200) 


27+ 1 ' 63 
"■^'-0.27 


3- u8 -3.08 





13 


0.33 



a Thc spectrum requires additional components for rcdshiftcd absorption edges originating from O VII or O VIII. 
b Thc inclusion of a Gaussian for the emission line at 6.4 kcV is required. 



Note. — (f) denotes a fixed value. Col.(l): Galaxy name. Col. (2): Observation start date. Col. (3): Photon index of the direct and indirect power-law components. Col. (4): Hydrogen column 
density within the source in units of 10 22 cm -2 . Col. (5): The power law dependence of the cmissivity (r~ e ). Col. (6): Inner radius of the accretion disk in gravitational radii with the outer 
radius fixed at the default value (100 gravitational radii). Col. (7): Inclination of the accretion disk to the line of sight in degrees. Col. (8): Iron abundance of the reflector relative to the Solar 
value. Col. (9): Ionization parameter of the gas. Col. (10): Reduced x 2 °f the best-fit model followed by the number of degrees of freedom in (). Col. (11): Observed 0.5—2 keV flux in units of 



10 12 erg s 1 cm 2 . Col. (12): Observed 2-10 keV flux in units of 10 12 erg s 1 cm 2 . Col. (13): Nominal absorption corrected 0.5-2 keV luminosity in units of 10 44 erg s 1 . Col. (14): Nominal 
absorption corrected 2-10 keV luminosity in units of 10 44 erg s _1 . 
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Table 6. PG QSOs: Constrained C06 Model a Best-fit Values 



Galaxy 


Obs 


. Date 




r 


N H 




i 


« 


xl (d.o.f.) 


Fo.5-2 


^2-10 


L 


5-2 


L-2 


- 10 






V z ) 




WJ 




y°) 






(8) 


( y ) 


( iu ) 


(11) 


(12) 


PG0050+124 


2002 


Jun 22 b 


2 


26 +0.02 


n 04+0-002 


4.5(f) 


30.0(f) 


10000(f) 


1.16 


(633) 


„ ,7+0.34 
23 


„ ,,7+0.24 
13 





98 





75 




2005 


Jul 18 b 


2 


16 +o.'oi 


i.9±g;l 


87.7^06 


28701^° 


1.05 


(810) 


5.16+° : fl 


4.87+iH 





55 





44 


PG0838+770 


2009 


Mar 2 


1 


4g + 0.'09 
— 0.09 




4.5(f) 


30.0(f) 


30(f) 


1.06 


(137) 


— 0. 03 


— 0.05 





14 





35 


PG0844+349 


2000 


Nov 5 


2 


32 +0.04 




1.7+0 6 


30.0(f) 


999+ 45 


1. 15 


267) 


6.76+3'^ 


4.5S+IH 





76 





45 


PG0953+414 


2001 


Nov 22 


2 


30+0.02 
— 0.02 




1 9+°'^ 
— 0.7 


30.0(f) 


48+jf 






3.39^°? 


2.93t° f 5 










PG1001+054 


2003 


May 4 




2.0(f) 




1.2 


30.0(f) 


297 


4.80 


(21) 














PG1004+130 


2003 


May 4 


1 


40 + 0.05 




4.5(f) 


30.0(f) 


30 +299 






o.io+°° 4 2 


0-33+°- 3 4 9 












2005 


Jan 5 


1 


94 +0.09 
z *-0.10 


14+ 008 
"■^-0.07 


4.5(f) 


30.0(f) 


30(f) 


1.10 


(961 


09+ 01 
u - ua -0.02 


43+ - 04 


o 


19 


o 


74 


PG1116+215 


2001 


Dec 2 


2 


„„ + 0.04 
°°-0.0i 




1 8+°' 8 


30.0(f) 


qqq + W 
yyy -410 


1.16 


(256) 


4 43+ 1 - 91 


2 97+ 1 - 64 


4 


03 


2 


57 




2004 


Dec 17 b 


2 


, 7 +0.01 
— 0.01 




— 0.6 


30.0(f) 


402 4 ; 32 


1.29 


(663) 


3.69+°'?, 4 
— 0.34 


3.09+0^°, 
— 0.20 


3 


35 


2 


68 




2004 


Dec 19 


2 


^ . +0 04 
1 -0.05 




, ~ + 9 
1D -0.4 


30.0(f) 


1433 4 :™" 


1.13 


(239) 


. ~~ + l 90 
4 20 T 
^' 2u -2.59 


oc + l -71 


3 


82 


2 


81 


ppii nfi H/1 1 












A Kff\ 

4.5(t; 






















PG1211+143 


2001 
2004 


Jviri 15 
Jun 21 




1.99 
2.02 




1.3 
1.5 


30.0(f) 
30.0(f) 


298 
300 


W21 
2.67 


(494) 
(539) 














PG1229+204 


2005 


Jul 9 b ' c 


2 


ns+0.04 
uo — 04 




4.5(f) 


30.0(f) 


391 ±59 


1.21 


(403) 


, , c+0.23 
olo -0.39 


, , c+0.23 
Jlo -0.37 





33 





30 


PG1244+026 


2001 


Jun 17 


2 


„ + 0.04 
09 




1.4+0-3 


87.3+ g 


299+ i:44 


0.86 


(238) 


6-61 + 386 


2.9S+° 1 ° 1 3 





39 





16 


PG1307+085 


2002 


Jun 13 b 


1 


— 0.02 




4.5(f) 


30.0(f) 


31+ 2 


0.98 


(145) 


— 0.08 


i.98+°';i 

— 0.23 





55 


1 


29 


PG1309+355 


2002 


Jun 10 


1 


61 + 0.02 
D1 -0.02 




4.5(f) 


30.0(f) 


3„3+0.29 


1.53 


(267) 


41+ - 11 
u -* 1 -0.07 


75+0- 17 
U - ' J -0.19 





04 





72 


PG1351+640 


2004 


Jun 23 b 


2 


oo + 0.03 
^-0.02 




4.5(f) 


30.0(f) 


30+2 


1.11 


(283) 


n sq+0.06 

U ' OM -0.03 


60+ 06 
U.DU_ Q q 4 





12 





12 


PG1411+442 


2002 


Jul 10 




2.0(f) 




4.5(f) 


30.0(f) 


285 


4.02 


(143) 














PG1426+015 


2000 


Jul 28 


2 


o,+0.07 
J1 -0.07 




i s+ 2 - 9 

iS -0.6 


30.0(f) 


30+28 


0.99 


(145) 


o n1 +l7.77 
8 - ul -8.01 


_ „„+14.77 
6 ' 67 -6.67 


1 


63 


1 


23 


PG1435-067 


2009 


Feb 2 


2 


2Q +0.13 
^ u -0.13 




4.5(f) 


30.0(f) 


67 +iio 

D '-37 


1.17 


(58) 


2.551°- 


2 60+ ' 80 
z - ou _1.06 


1 


29 


1 


09 


PG1440+356 


2001 


Dec 23 


2 


76 + 0.04 
,D -0.04 




1 7+ 1 - 5 


30.0(f) 


30+20 


1.21 


(163) 


5.19™ 


2 21+ 4 ' 12 
— 2 21 





83 





34 




2003 


Jan 1 


2 


Ko + 0.02 
°°-0.02 




L - '-0.4 


30.0(f) 


30 +1 


1.67 


(434) 


r rr + 2.29 

o.oo_ 1 66 


9 o n +l.'26 
2-»°-0.90 





89 





43 




2003 


Jan 4 


2 


74 +0.02 
' -0.02 




1 6+°' 8 


30.0(f) 


30+3 


1.27 


(343) 




r* i o + 2 . 25 

2,1 -2.13 





79 





33 




2003 


Jan 7 




2.0(f) 




1.2 


30.0(f) 


622 


3.93 


(340) 














PG1448+273 


2003 


Feb 8 


2 


.,,+0.04 




, 7 +0.6 


30.0(f) 




1.37 


(382) 


3.97+lH 


1 81+ 1 - 13 
1 ' Si -0.83 





45 





19 


PG1501+106 


2001 
2001 


Jan 13 
Jan 13 


2 


,0+0.01 
^ a -0.01 
2.27 




1 8+°' 3 
1.7 


30.0(f) 
30.0(f) 


3 « +1 

958 


1.62 
2.57 


(484) 
(570) 


14 - 23 ±0.29 


12.62+ ; 7 3 





47 





38 




2001 


Jan 14 


2 


, fi +0.01 
' !D -0.01 




l s+ - 3 
18 -0.6 


30.0(f) 


3 « +1 


1.44 


(627) 


i6.i8t;;S? 


14 -82t ;^ 





53 





45 




2005 


Jan 16 




1.51 




1.7 


30.0(f) 


289 


3.96 


(825) 
















2005 


Jul 17 b 


1 


63 +0.02 
DJ -0.04 




1 6+ - 4 


30.0(f) 


488+ 88 


1.53 


(693) 


K «« + 0. 65 
6-88_ g7 


11 46+ 1 - 10 





57 





35 


PG1613+658 


2001 


Apr 13 


1 


qq+0.14 
3 -0.14 


23+ 10 
u ^ J -0.20 


j 2 + 53.6 


50.6+1.1. 


no+ 7 f 


0.97 


(132) 


2.49+1^ 


4 19+ 5 - 64 
^■- L -2.10 


2 


22 


1 


87 




2001 


Aug 29 


1 


„„ + 0.15 
— 0.08 


o oi +0 - 13 

u,ul -0.01 


4.5(f) 


47 3+ 7 '° 


30+1089 


1.00 


(244) 


4.08±itS 


5 26+ 1 ' 41 
J ' ZD -0.90 


1 


97 


2 


31 


PG2130+099 


2003 


May 16 




1.59 




1.6 


30.0(f) 


300 


2.86 


(485) 


2-5^0.72 


4.65^;^ 










B2 2201+31A 


2008 


Dec 1 


2 






2.2 4 J 4 1 ° 


30.0(f) 


sit; 


0.97 


(213) 


10.08 


13 


.02 


PG2214+139 


2002 


May 12 




2.0(f) 




4.5(f) 


30.0(f) 


30 


8.42 


(210) 















a Thc accretion disk emissivity index of the kdblur component and the iron abundance of the rcfiion component arc fixed at 3.0 and 
Solar, respectively. The error estimates of individual parameters, fluxes, and luminosities for poor fits (i.e. xf y > 2.0) arc not calculated. 

b The spectrum requires additional components for redshifted absorption edges originating from O VII or O VIII. 

c Thc inclusion of a Gaussian for the emission line at 6.4 keV is required. 



Note. — (f) denotes a fixed value. Col.(l): Galaxy name. Col. (2): Observation start date. Col. (3): Photon index of the direct and 
indirect power-law components. Col. (4): Hydrogen column density within the source in units of 10 22 cm~ 2 . Col. (5): I nncr radius of the 
accretion disk in gravitational radii with the outer radius fixed at the default value (100 gravitational radii). Col. (6): Inclination of the 
accretion disk to the line of sight in degrees. Col. (7): Ionization parameter of the gas. Col. (8): Reduced x 2 °f t- ne best-fit model followed 
by the number of degrees of freedom in (). Col. (9): Observed 0.5-2 keV flux in units of 10~ 12 erg s~ 1 cm -2 . Col. (10): Observed 2-10 keV 
flux in units of 10 — 12 erg s — 1 cm -2 . Col. (11): Nominal absorption corrected 0.5—2 keV luminosity in units of 10 44 erg s — 1 . Col. (12): 
Nominal absorption corrected 2—10 keV luminosity in units of 10 44 erg s — 



-40- 



Table 7. Hardness Ratio Estimates of Weak''" ULIRGs 



Galaxy 


Date 


Extended? 


Total 


Hard 


Soft 


HR 


r 


N H 


^0.5-2 keV ^2-10 keV L Q.5-2 keV L 2-10 keV 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) (12) (13) 



F00091 — 0738 


2008 


Nov 1 


No 


7 


, + 3. 

3-1 


6 
6 




2 
.9 


-0 


14+0-59 
- L -0.33 


1.03+; 


66 
.13 


„ , , + 1.47 
°-44_ .40 





1 1 


0. 


56 





.04 





.20 


F00188-0856 


2003 


Sep 4 


No 


16 


„-|-3 

D -2 


6 
.4 


10+4 


3 
.1 


-0 


, r + 0.36 
-0.25 


l.llu 


5 
6 


,-. „,>4-0 64 

o.36lg;°t 





.42 


1. 


69 





20 





.71 




2004 


Dec 20 




164 


,-.-,-4-8 
52 _7 


8 
.7 


112 +10.7 
-9.7 


-0 


, 7+0. 00 
-0.08 


1.45+° 


16 
.17 


,4-0 09 
UJ ^-0.02 





.42 


1. 


12 





.20 





.47 


F00456— 2904: SW 


2009 


May 22 


Yes 


43 


,,-1-4 

y -2 


1 

.9 


34+° 


9 

.8 


-0 


91 +0.22 
" !i -0.17 


1 15+° 
iio_ 


.31 
.41 


^ ^~-r-0 35 
% 5 ' ' 





.74 


3. 


11 





.22 





.94 


F01004-2237 


2003 


Aug 3 


No 


20 


6I2 


6 
4 


14+ 4 
1 -3 


8 
.7 


-0 


4Q +0.32 
' w -0.24 


1 4+° 


6 
6 


0.16+ ' 47 





.26 


3. 


64 





.09 


1 


.29 


F02021-2103 


2009 


Oct 12 


No 


39 


- ,-,-1-5 

18+° 


3 
2 


21+.1 


.7 
.6 


-0 


„„ + 0.20 
08 -0.16 


0.89+g 


.29 
.36 


^^4-0 36 
U - O0 -0.23 





56 


3. 


56 





.19 


1 


21 


IRAS03521+0028 


2002 


Dec 25 


No 


3 


- 4-2 


3 

.8 


2+_\ 


.6 
.3 


-0 


,, + 1.23 
•"-0.54 


I.5+3 


9 
.4 


0.25+ 5 ' 64 





14 


0. 


33 





.09 





20 


F04103-2838 


2003 


Apr 28 


No 


30 


12+I 


6 

.4 


18l 5 4 


.3 
2 


-0 


2 „ + 0.24 
zu -0.18 


1.06±g 


35 
.45 


^ . ^4-0 40 
(,m ^°-0.24 





.74 


3. 


37 





.26 


1 


.17 




2006 


Feb 13 


No 


184 


79 -8 


9 
9 




3 
.2 


-0 


14 +0.08 
-0.07 


i.oelg 


11 
12 







.80 


4. 


10 





.28 


1 


.42 


F05024-1941 


2007 


Feb 8 




129 


573°' 




72+_% 


.5 
.5 


-0 


,,+0.10 
-0.09 


1 05+° 
i.uo_ 


.14 
.17 


°- 40 ±0.12 





.79 


3. 


91 





81 


4 


.02 


F07598+6508: E 


2001 


Oct 25 




240 


En +8. 
50 _7 


1 



190±£ 


9 
.8 


-0 


7n +0.08 
'"-0.07 


2.261? 


24 
.23 




2 


.45 


1. 


93 


1 


.42 


1 


.12 


F08572+3915: NW 


2006 


Jan 26 


No 


9 


-7+3 
' -2 


§ 
6 


9+2 
^-1 


.6 
.3 


0.56+°;« 


-0-43±°J 7 







.06 


3 


01 





.01 





24 


F09039+0503 


2008 


Dec 31 


Yes 


34 


-4-3 
' -2 


8 
6 


27+» 


.3 
2 


-0 


59 +0.25 
J -0.20 


2.00+g 


.70 
.57 


<0.18 





.58 


0. 


68 





23 





27 


UGC 5101 


2001 


May 28 


Yes 


512 


, « .4-15 




9 


318+i? 


9 
.8 


-0 


24 +0.05 
z -0.05 


i.iolg 


10 

.09 


„ m „ + 09 

o.selgg? 


1 


61 


7. 


36 





.06 





25 




2001 


Nov 12 




698 


318+^' 


9 
8 


380+ 2 g 


5 

.5 


-0 


09 +0.04 
u -0.04 


0.97+g 


07 
.07 


^ . ^4-0 05 
u - 4O -0.07 


2 


.67 


14.93 





.09 





.51 


F10190+1322 


2003 


Jan 31 


No 


16 


6+1 


6 
4 


10+ 4 


.3 
.1 


-0 


,5 + 0.36 
ZO -0.25 


1.18lg 


50 
.68 


n oe + 0-65 
u.oo_ Q 3Q 





.46 


1 


70 





.07 





.24 




2003 


May 5 




131 


45+1 


8 
.7 




.3 
i 


-0 


,.+0.10 
Ji -0.09 


l.S6±° 


.19 
.16 


n 1Q +o.io 

u y -0.08 





97 


2. 


99 





. 14 





42 


F10565+2448 


2003 


Oct 23 


Yes 


346 


»iti°, 


.6 


,55 + 17.0 
^ 00 -16.0 


-0 


47 +0.06 
*'-0.06 


1.62+° 


14 

.13 


^ ^,^4-0 07 
°- 05 l0.04 


2 


.43 


5. 


00 





.10 





.21 




2003 


Jun 17 




691 


166m 


9 
9 


525+ 22 


9 
.9 


-0 


c-9 + 0.05 
' D -0.04 


i.nlg 


09 
.11 


<0.56 


3 


.83 


6 


87 





.16 





.29 


F11095-0238 


2009 


Apr 9 


No 


48 





9 


33± 6 5 


.8 
.7 


-0 


,o + 0.19 
JS -0.15 


1.63±g 


.34 
.39 


^ - „+n 24 
0.13 + 





.77 


1. 


81 





22 





.52 


F12072 — 0444 


2003 


Feb 1 


No 


16 


2 +2. 


6 
3 


"±S 


8 
.7 


-0 


7^ + 0.43 
,0 -0.25 


2.5+? 


6 
1 







.66 


0. 


36 





28 





15 


F12112+0305 


2001 


Dec 30 




227 


77+_l 


8 
.8 


150i« 


3 
.3 


-0 


, + 0. 08 
' oZ -0.07 


1.33+g 


12 

.14 


18+2^ 

u c -0.06 


1 


.34 


4. 


31 





.17 





.55 




2003 


Apr 15 


No 


51 


14+t 


8 
.7 


37+ 7 


.1 
.1 


-0 


45 + 0.19 
*°-0.15 




4 
.4 


0.11+ ' 24 


1 


50 


3. 


26 





19 





.41 


F13218+0552 


2004 


Jul 11 




41 


15+ 5 ' 



9 


26+f 


2 
.1 


-0 


, 7 +0.20 
■"-0.16 


1.39lg 


.29 
.33 


17+ - 21 
' -0.05 


1 


29 


3. 


77 


1 


.54 


4 


.49 


F14348-1447 


2002 


Jul 29 




117 


31+= 


6 


86±i°, 


.3 
1 


-0 


,7+0.12 
' -0.10 


1 77+° 
L - " -0 


22 
.24 


<0.15 


2 


73 


4. 


47 





.45 





.74 


F14348-1447: NE 


2006 


Mar 12 


Yes 


24 


6+1 


6 
4 


18tj' 


.3 
.2 


-0 


+ 0.25 
ou -0.17 


1.83lg 


.62 
.69 


<0.32 





.78 


1. 


17 





13 





19 


F14348-1447: SW 


2006 


Mar 12 


No 


33 


13+4 


7 

6 


20±5 


.6 
.4 


-0 


,,+0.22 
zl -0.17 


1.20lg 


.43 
.52 


37+ - 37 
Uii '-0.21 





.77 


3. 


02 





.13 





.50 


F15130-1958 


2003 


Jun 2 


No 


38 


7+1 


8 
6 


31+| 


.6 


-0 


6 , + 0.24 
Do -0.19 


2.15lg 


.75 
.65 


<0.17 


1 


.59 


1. 


33 





.47 





40 


F15250+3608 


2002 


Feb 22 




215 


38+1 


2 
.2 


177+ 14 
11 '-13 


3 
3 


-0 


„ 5 + 0.09 
° -0.08 


2.02lg 


24 

.23 




1 


.76 


1. 


98 





.12 





.14 




2003 


Aug 27 


No 


37 


5 +3. 
°-2 


4 
.2 


32l 6 5 


.7 

.6 


-0 


7, + 0.25 
' -0.20 


2.271* 


24 
.77 


<0.07 


1 


.43 


0. 


96 





10 





07 


F15462-0450 


2009 


Apr 23 


No 


496 


260tiJ 


1 
1 


236+ 16 4 
^ JD -15.4 


05+ 05 
u ' uo -0.04 


0.78+g 


09 
.09 


79+°- 10 
u -' y -0.09 


7 


.02 


52.70 


1 


74 


13.05 


F16090-0139 


2003 


Feb 10 


Yes 


27 


10+ 4 ' 


3 
1 




2 
.1 


-0 


+ 0.23 
-0.20 


1.57lg 


53 
.45 


0.15+ 033 





.82 


2. 


34 





.38 


1 


.09 


F17208-0014 


2001 


Oct 25 


Yes 


478 


155+_\l 


5 
4 


323li| 



.0 


-0 


,5 + 0.05 
JO -0.05 


1.49lg 


.10 
.11 


21+ 06 
u -' al -0.06 


2 


.12 


5. 


31 





.09 





.22 




2002 


Feb 19 




245 


o,+ 10.1 


1631 1 ! 


8 
.8 


-0 


,, + 0.07 

•"-0.07 


i.sslg 


.13 
.13 


17+ - 06 


3 


.18 


7. 


29 





.13 





.31 




2003 


May 7 


Yes 


92 


30l 6 5 


5 

.5 


62+ 8 7 


9 
.9 


-0 


,5 + 0.13 
' ao -0.11 


1.43lg 


27 

.23 


23+0.18 


3 


.52 


8. 


49 





15 





.36 


F21208-0519: N 


2009 


Mar 31 


No 


14 


, + 2. 


9 
6 


11+ 4 
1 -3 


4 

.3 


-0 


57+0.43 
° ' -0.30 


2.001J 


.20 
.95 


<0.45 





.33 


0. 


38 





.14 





.17 


F21208-0519: S 


2009 


Mar 31 


No 


6 


1 + 2 


3 

.8 




4 
.2 


-0 


fi7 +0.82 
D '-0.47 


23li 


8 


<1.00 





.17 


0. 


12 





.07 





.05 


F21329-2346 


2009 


Jun 21 


No 


27 


9+J- 


1 

9 


18l 5 4 


.3 
2 


-0 


,, + 0.26 
aa -0.20 


1 40+ 44 
1 -* u -0.50 


20+ ' 36 





51 


1 


47 





20 





.59 


F22491-1808 


2001 


May 24 




174 


34+6 


9 

.8 


140+1= 


9 
.8 


-0 


61 +o.io 

,DJ --0.09 


1.96+g 


2 b 
.23 




1 


16 


1. 


43 





.17 





.21 




2007 


Jul 13 


Yes 


41 


5 +3. 

°-2 


4 

2 


*»±l 


1 

.0 


-0 


7g+0.15 
' ,D -0.10 


2.50lg 


60 
.61 




1 


.13 


0. 


62 





.17 





.09 


F23234+0946 


2009 


Aug 15 


Yes 


39 


x + 2. 
-0 


3 
8 


38l^ 


2 
.1 


-0 


,5 + 0.27 
ao -0.05 


4.0lj 


.0 




1 


.09 


0. 


08 





.46 





03 



Note. — t We arbitrarily consider a source "weak" when its count rate is < 0.01 and 0.05 counts per second when observed with Chandra and XMM-Newton, 
respectively. Col.(l): Galaxy name. Coordinate-based names beginning with "F" are sources in the IRAS Faint Source Catalog. Col. (2): Observation start date. 
Col. (3): Whether the source is extended when compared with the telescope PSF. The symbol " ■ • • " denotes XMM-Newton observations where radial profiles 
were not created; all these observations are assumed to be unresolved due to the large PSF of the telescope relative to Chandra. Col. (4): Total counts in the 
0.5-8.0 keV band. Col. (5): Total counts in the 0.5-2.0 keV band. Col. (6): Total counts in the 2.0-8.0 keV band. Col. (7): The hardness ratio. Col. (8): The 
estimated T based on the hardness ratio and assuming Galactic absorption. Col. (9): The estimated total absorption assuming T — 1.7 in units of 10 22 cm 2 . 
Col. (10): The estimated observed 0.5—2.0 keV flux assuming the estimated T and Galactic absorption in units of 10 — 14 erg s 1 cm 2 . Col. (11): The estimated 
observed 2.0—10.0 keV flux assuming the estimated F and Galactic absorption in units of 10 14 erg s 1 cm 2 . Col. (12): Derived 0.5-2 keV luminosity from 
Col. (10) in units of 10 42 erg s _1 . Col. (13): Derived 2-10 keV luminosity from Col. (11) in units of 10 42 erg s _1 . 



Table 8. Best-fit Parameters for ULIRG Spectra with > 100 Counts 



Galaxy Date Model kT T N n Bource i / COV erl N n Bourco2 /covor2 E, inr a EW Stat 



-lOkcV 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) 



Cash Statistics (unbinned spectra with < 1000 counts) 



F00188-0856 


2004 


Dec 20 


A 








1 


76+° 
— 


33 
29 


















716.0 


1940 





42+g 
— 


16 
18 


0.77+gS 
— 0.63 


0.20 


0.33 


F04103-2838 a 2006 


Feb 13 


C 





10+g 


03 
08 


1 


38_° 


97 

;H 


19+°- 33 
uiB -0.19 












1.62+1 


58 
12 


848.1 


1953 


2 


17+° 
1 -0 


87 
87 


4 56+ 2 ' 04 
° -0.97 


0.64 


1.35 


F05024-1941 


2007 


Feb 8 


D 








1 




47 


32.9+29 3 5 9 





86+.° 


40 










646.5 


1938 







50 


o o« + 2 . 10 
2 ' 2B -2.28 


4.45 


5.62 


F07598+6S08 


2001 


Oct 25 


D 








2 




27 
23 


„ ,, + 44.8 





95+° 


05 
15 










655.6 


1938 


2 


02+2 


61 
02 


2 50+ 2 98 


29.17 


8.40 


F10190+1322 


2003 


May 5 


A 








1 


78+° 


62 
31 


















661.5 


1940 





62+.°, 


30 
34 


1 11 + 12 ° 
i,ii -1.18 


0.10 


0.16 


F10565+2448 


2003 


Oct 23 


C 





68i° 


17 
12 


2 


14+° 
1 -0 


66 
53 


30+ 57 

U.3U_ Q 21 
















352.3 


507 


3 


19+ 2 


00 
06 


, 7o + 2.07 
J ' '^-4.18 


0.23 


0.17 




2003 


.Tun 17 


C 





68±° 


14 
07 


1 


41+° 


25 
23 


03+°'° 4 
















1118.9 


1937 


3 


«±°0 


67 
73 




0.17 


0.37 


F12112+0305 


2001 


Dec 30 


B 





11+g 


07 
03 
























581.9 


1940 





l°t°o 


20 
10 


0.00+ 119 


0.01 


0.00 


F14348-1447 


2002 


Jul 29 


A 








1 


97+g 


45 
29 


















521.2 


1940 


1 


62+.°, 


78 

82 


2.48+ 2 2 3 2 ° 


0.35 


0.43 


F15250+3608 


2002 


Feb 22 


D b 







12 
18 


2 




85 
43 


, ,, „ + 202.0 
















613.2 


1936 


1 


54t°o 


42 
48 


3.54ti:°3 2 8 


11.31 


5.65 


Arp 220 


2000 


Jun 24 


C 







05 
05 


1 


82+° 


19 
19 


125.2+_|; 2 





98+° 


01 
04 


1.83(f) 


0.0(f) 


0.05+° 


22 
05 


605.4 


505 


6 


31+1 


59 
99 




1.69 


2.55 




2002 


Aug 11 


c 





°*±°0 


12 
07 


1 


65t°o 


17 
15 










6.70(f) 


0.0(f) 


1.2115 


40 
21 


1213.6 


1937 


7 


24+5 


09 
04 


11 58+ 350 
11.00_ 3 46 


0.06 


0.08 




2003 


Jan 15 


c 





63+.U 


08 
18 


1 




19 
23 










6.70(f) 


0.0(f) 


o.8i_J 


20 
81 


1019.7 


1937 


6 


51±i 


16 
23 


91 18+ 4 ' 57 
zi ' io -10.03 


0.05 


0.12 




2005 


Jan 14 


c 





62_° 


08 
12 


1 


«tg 


39 
50 


7 o ,+136.5 
' al -63.0 





88t° Q 


10 

36 










849.7 


1936 


6 


16+3 


94 
46 


, - 7^ + 8. 13 
10 -' D _15.74 


0.22 


0.50 




2005 


Feb 20 


c 





5 y ±°o 


11 

24 


1 


"t°0 


54 
36 


















585.7 


1248 


7 


72+. 2 


21 
26 


17 03 + 10 89 
LIU3 -17 .03 


0.06 


0.12 


F15462-0450 


2009 


Apr 23 


D 




74t° 


14 
16 


1 


»t°0 


38 
37 


1 30+°- 78 





79+° 


10 
20 










479.6 


508 


6 


i»t; 


42 
19 


41 06+ 9 ' 55 


4.92 


1 — 1 

11.17 


F17208-0014 


2001 


Oct 25 


C 





2 


44l 


63 
88 


43+ 83 
u '^ a -0.30 








1.83(f) 


0.0(f) 


07+" a5 
u ' u '-0.07 


484.5 


506 


1 


53±°. 


77 
40 


1 59+ 1 ' 73 
loy -1.59 


0.16 


0.07 




2002 


Feb 19 


C 







16 
13 





°0±°0 


56 
53 


















863.1 


1938 


1 


o, + 


68 
73 


7 98+ 4 53 
'' yo -7.63 


0.11 


0.35 


F22491-1808 


2001 


May 24 


C 





7lt°o 


17 
16 


2 


20+° 


42 
38 


















633.4 


1938 


1 


17+° 
i '-0 


35 
33 


69 + 123 
U - Dy _0.69 


0.19 


0.10 



X 2 Statistics (spectra binned to at least 15 counts bin 1 ) 



E 


q KK +1.02 
2-66_ 66 


2.54+_° 


09 
09 










C 


3+ ' 1 


2 .1±° 


1 
1 










D 


°- 08 -o.oi 


2.23t° 


06 
07 


7 




53 
51 


n qq+0 -01 
°- 98 -o.oo 


D 


20+°' 04 

U -2U_ Q Q2 


lA9+° 


10 
22 


6 


67lg 


70 
75 


96+°' 3 ° 
u - yD -0.33 


D 


„ 21 + 0.04 


1.41+g 


14 
22 


5 


8lt8 


68 
88 


97+°- 01 
u,y ' -0.01 


D 


0.70+g;» 


1.27+° 


58 
35 





l 7 t°o 


50 
15 


1.00(f) 


D 


o.85±° ; 3 ? 


l.8l+° 


40 
40 





1»±S 


16 
09 


1.00(f) 


D b 




1.97+° 


17 
16 


1 


71+° 


20 
19 




E' 1 


0.94+™; 


2.07+_° Q 


01 
01 











PG0157+001 2000 Jul 29 

2003 Aug 23 
F05189-2524 2001 Mar 17 

2001 Oct 30 

2002 Jan 30 
UGC 5101 2001 May 28 

2001 Nov 12 
F11119+3257 2002 Jun 30 
3C 273 2000 Jun 14 

2000 Jun 17 

2001 Jun 13 

2003 Jul 7 

2004 Jun 30 

2007 Jan 12 

Mrk 231 2000 Oct 19 
2001 Jun 7 
2003 Feb 3 



137.66+*?;*! 0.96+8;% 

78.23± 32 ; 32 7 o.97±° ;° 2 



2 94+°' 15 
z - y -0.14 

0.08+ ' 002 



-0.002 
s+0.11 
-0.10 
q+0.02 
-0.01 
,+0.05 
-0.04 
q 7Q+0.04 
u - /y -0.03 

2 91+°- 13 
" iyi -0.11 
„ OO+0.07 
°- 82 -0.06 
,+0.24 
-0.18 

93+°-° 4 
u - y,3 -0.04 
, ,,+0.01 

°- i "-o.oi 

69+°' 2 ° 
U - Dy -0.11 
c+0.11 
-0.08 
n+0.10 
-0.09 



2.652 
0.891 
2.93 4 



2. 84"' 



0.49: 
0.80 4 



,+0.01 
-0.01 



9 9S +0.01 
2 - 28 -0.01 



2 08+° 01 
z - uo -o.oi 



„ .,+0.02 
z - 1,3 -0.02 



C + 0.02 
°-0.01 



1.44 
1.87 



+ 
-0 


46 
23 


3 


°3t. 


61 
99 





7 g+0.12 
'°-0.17 


+ 
-0 


07 
04 


8 


3415 


89 
71 





7R +0.03 
,D -0.04 


+ 
-0 


40 
20 


3 


75^5 


09 
17 





„ n +0.10 
S0 -0.11 



7 +25.77 
-12.45 



n+o.io 

"-0.12 



136.4 
138.2 
164.1 
131.8 
114.6 
28.6 
38.2 
97.1 



148 
130 
165 
118 
108 
24 
45 
106 



2.28 
21 

1587.5 1571 4.30 



101 69 + 8 - 81 

rn „ + 14.24 
150.16_ 25 n6 

I2.i9l 2 ; 9 5 
9.46l 2 ;° 6 2 

9 49+ 1 - 67 
y y — 4.30 

2 04+ 4 95 
z,u -0.74 
+ 2.91 
— 1 99 

26.411 4 ; 29 

+0.02 
-0.02 



X 10 3 8.08 



, + 27.06 
-42.53 
K+12.14 
-11.96 
302.82l^; 7 8 

256.44lf 09 9 44 

336.9ll 4 ? 7 49 8 

16.781? 7 ;? 7 

I9.45l5 3 9 ; 4 ° 3 

155 91+ 17 49 
+0.15 

0.14 



96.742 
134. 76^ 



X 10 3 



1600.6 1481 1.00+° °5 X 10 3 7.72+" 09 x 10 3 2.89X10 3 



1922.4 1725 5.39+g'g 2 X 10 3 9.66+°'5o X 10 3 3.90X10 3 



1701.2 1589 6.38+° ° 3 X 10 3 9.79+°'5 3 X 10 3 4.60X10 3 



978.7 943 4.0ll° ° 3 X 10 3 6.80+g'?- 1 , X 10 3 2.89X10 3 



1843.7 1708 4.87 



+0.03 
0.03 



7 +0.03 
'-0.02 



+ 0.14 
-0.13 



.+0.14 
^-0.10 



+ 1.15 
"-0.40 



87.7 
139.4 
82.1 



.74"' 



133 
91 



*t 7 : 2 67 

I0.22l 3 ;°? 
9.08+i' 62 



R o , 7 +26.19 
6S.17_ 67 14 

n + 19.07 

-33.93 



63. 50^ 



61 54+ 16 25 



83.20 


67.00 


99.99 


72.99 


28.06 


23.04 


6.41 


15.85 


7.23 


19.96 


1.85 


6.45 


3.06 


4.75 


143.92 


173.32 


3.10X10 3 


5.52 xl() 3 


2.89 X10 3 


5.27X10 3 


3.90X10 3 


6.60X10 3 


4.60X10 3 


6.69X10 3 


2.89X10 3 


4.65 xll) 3 


3.50 X10 3 


8.73X10 3 


4.95 


7.88 


1.39 


3.54 


4.25 


5.82 



Table 8 — Continued 



Galaxy 


Date 


Model 


kT 


r 


JV H , soured /covorl A r H, B ourcc2 /covor2 




a 


EW 


Stat d.o.f F 5_2kcV 


^2-lOksV L 0°5-2kcV L 2-TokcV 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) (7) (8) (9) 


(10) 


(ii) 


(12) 


(13) (14) (15) 


(16) (17) (18) 





2003 


Feb 


11 


D 


0.65 




2003 


Feb 


20 


B 


0.66 


Mrk 273 


2000 


Apr 


19 


D 


0.86 




2002 


May 7 


D 


0.70 


F13451+1232 


2000 


Feb 


24 


D 




NGC 6240 


2000 


Sep 


22 


D 


0.76 




2001 


Jul 


29 


D 


0.86 



+ 0.05 
-0.05 
+ 0.05 
-0.05 
+ 0.11 
-0.07 
+ 0.07 
-0.05 



1 16 + 019 
25+ 07 
1.63 



o + 3.14 
-2.32 



+ 0.16 
0.18 



41. 29 1 



+ 6-03 

.86+^ 59.32±*1 
+ 0.12 
0.14 



1.55"' 



2.14 



+ 0.13 
-0.12 



11.54 
, Q . + 0.53 

19+0-03 
u -0.03 



. + 0.12 
3 -0.14 



97+ 01 
u ' a '-0.01 

1.00(f) 



K +14.32 
D — 12.71 



,+0.01 
-0.01 



+ 0.03 
-0.03 



3 + 0.15 
J -0.13 



.+0.04 
-0.04 



1.00(f) 



bb.49_ 2g 56 0.81_ 06 



2002 Mar 12 D 0.62+°'™ 2.04+° ^ 0.20+°\°l 1.00(f) 142.69tJ 7 ; 2 3 °' 98 -0.01 



2003 Mar 14 D 0.821"-"' 2.001" *-' 0.21 



+0.07 2 Qo+0.25 „ 21 +0.07 1Mffl 109 25+ 34 " 12 94+ 03 
-0.06 -0.21 u - z -0.06 t.uu^i; iuy.zo_ 22 24 u.y^_ Q 05 



2003 Aug 13 



D o.82±S; ;« 1.82+2:3 0.14±° °$ 1.00(f) 146.«±*;£ 0.96±° ; ™ 



2003 Aug 21 D 0.82^009 l- 90 t u .20 °- 16 -o'.05 l-00(f) 201.58jl||;U 0.99+.° 01 



2003 Aug 29 



, + 0.15 , .,+0.32 „ n7 +0.09 



.,+48.72 n Qfi + 0.01 



+ 0.03 
-0.03 
+ 0.15 
-0.16 



1.83(f 
6.40(f 
6.67(f 
7.01(f 
1.83(f 
1.94(f 
2.40(f 
6.40(f 
6.67(f 
7.01(f 
1.06(f 
1.83(f 
6.40(f 
6.67(f 
1.83(f 
6.40(f 
6.67(f 
7.01(f 
1.83(f 
6.40(f 
6.67(f 
7.01(f 
1.83(f 
6.40(f 
6.67(f 
1.38(f 
1.83(f 
6.40(f 



0.00 
0.25 



+ 0.06 
0.00 
+ 0.20 
0.11 



0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 
0.00(f 



0.23 
0.67 



+ 0.62 
0.00 
+ 4.44 
0.67 



06+ 05 
u - UD -0.04 
n or+0.18 

n+0- 12 
u,11 -o.n 

o.oo+°- 20 



. 00 
+ 0.24 
0.04 
+ 0.53 
0.03 
+ 0.91 
0.06 
+ 0.47 
0.46 
0.31 
0.14 
+ 0.85 
0.13 
+ 0.05 
0.03 
+ 0.07 
0.04 



0.04 
0.07 
0.09 
0.46 
0.20 
0.13 
0.05 
0.04 

17+0-23 
' -0.17 

09+ - 20 
u - ua -0.09 

n ns+0- 42 
0-08-0.08 

„ , 9 +0.39 
u - zz -0.22 

17+0- 39 
-0-17 

09+ - 96 
u - UM -0.09 

04+ 05 
u - u -0.04 

25+ - 21 
u -^°-0.25 

14+0- 20 
u - 1 *-0-14 

0.03+ - 24 



0.033 

0.52j 
0.07 4 



0.123 

3 +l 

-0.20 



101.3 
76.31 

67.3 
115.4 

53.4 
365.8 



95 
7!) 
85 
89 
60 
333 



B-15i°l 9 07 
7.08 



+ 0.76 
0.77 
1.00 



9.07 
62.64 



+ 2.70 
5.98 



66 , e +23. 75 
on ,,+ 10.34 
82.12+.^;™ 
37 - 2 °l22".?5 

io 2 -«±72-.™ 

188.90± 44 ™ 



250.9 203 54.98 



+3.55 lfi , 0.+21.39 
36.31 ± " ao -147.68 



171 56.02+. 4 - 20 , 199.13± 4 ^; 4 | 



102 55.38±I£i 4 190.13±« 4 7 47 4 



264.1 237 56.34 



-0.03 
, + 0.19 
-0.03 
10.38 
-0.52 
7+0.47 
-0.07 

5 + 1-14 
-0.12 
0.20+^ 6 n 7 ' 42 



+ 120-16 o ne - HR + 269.44 



111 56 91~^ 5 ' 30 211 qi+ 25 - 69 



51 58.15™/ 227.6S; 342 - 2 



1.00 
0.35 
5.28 
4.29 
21.56 
60.70 



3.56 
3.31 
10.73 
6.44 
49.05 
59.78 



+- 

to 



"Values taken from iTeng et aTl )2008l) . The energy bands in that paper are 0.2—2 keV and 2—10 kcV. 

^The ratio of the normalization parameters between the scattered and direct components is 249 for F15250 + 3508, 62 for F1119+3257, 24 and IS for the 2000 and 2002 observation of Mrk 273, respectively. 

c Thc addition of an emission line at 6.4 kcV docs improve the fit, but not statistically significant. 

d The addition of a red-shifted oxygen edge improves the fit significantly. 

c Thc addition of this absorption component docs improve the fit, but not significantly. 

Note. — Col.(l): Galaxy name. Coordinate-based names beginning with "F" are sources in the IRAS Faint Source Catalog. Col. (2): Observation start date. Col. (3): Best-fit continuum model, all 
modified by Galactic absorption: A — single power law, B — MEKAL only, C — MEKAL plus power law, D — scattering model or partial covering model, E — blackbody plus power law. Col. (4): MEKAL 
temperature, if applicable, in units of kcV. Col. (5): Spectral index of the power law (or reflection) model. Col. (6)— (9): Absorption within the source and corresponding covering fraction for the covering 
fraction absorber model. The absorption values are in units of 10 22 cm -2 . Col. (10): Rest energy of the emission line(s) in units of keV. Col. (11): Line width(s) of the emssion line(s) in units of keV. 
Col. (12): Equivalent width(s) of the emission linc(s) in units of keV. Col. (13)— (14): Fitting statiscs and degrees of freedom, respectively. Col. (15)— (16): Observed flux after background subtraction in 
units of 10 — 14 erg s — 1 cm" -2 . Col. (17)— (18): Nominal absorption corrected luminosities in units of 10 42 erg s — . 
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Table 9. PG QSOs: Tested Universal Spectral Models 



ID 


AGN 


SE 


xl 


(1) 


(2) 


(3) 


(4) 


A 


PL 




8.0 


B 


PL + AbsxPL 




3.0 


C 


PCFxPL 




3.4 


D 


PCFxPCFxPL 




1.8 


E 


PCFxPCFxPL 


MEKAL 


1.5 


F 


PL 


CompTT 


2.1 


G 


PL 


REFLION (blurred; C06) 


19.9 


H 


PL 


PEXRAV 


3.6 


I 


PL 


PEXRIV 


3.1 


J 


PL+PEXRAV 


zBbody 


1.3 


K 


PL+PEXRIV 


zBbody 


1.2 


L 


C06 


zBbody 


1.4 


M 


PL+PEXRIV 


zBbody +Zcdgc 


1.1 



Note. — Col.(l): Model identifier. Col.(2): Spec- 
tral models fitting the AGN component: PL = red- 
shifted power law, AbsxPL = absorbed redshiftcd 
power law, PCF = partial covering absorber. Col. (3): 
Spectral models fitting the soft excess (SE) com- 
ponent: CompTT = Comptonization model, RE- 
FLION= blurred ionized reflection model of C06, 
PEXRAV = neutral reflection, PEXRIV = ionized 
reflection, ZBbody = redshiftcd blackbody, Zedge = 
redshifted absorption edge. (4): Reduced \ 2 from the 
fits. 
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Table 10. PG QSOs: Derived Parameters from Model M 



Soft Excess Power Law Reflection 



Galaxy 




K 


^0.5— 2 


Lo.5—2 


^0.5-2 


^2 — 10 


Lo.5—2 


£2-10 


Fo.5-2 


^2 — 10 


^0.5 — 2 


-^2 — 10 








1 v i n- 14 


\ /yin 42- ! 

) (X1U ) 


l v i n— 12 \ 

(X1U ) 


^ X ±u 


) l xlu 


) (xlO 44 ) (xlO- 14 ) 


^ X 1U 


l v 1 n 44 


) l xiu 




o\ 
\ z ) 


l%\ 


IA\ 

\V 


(°) 


v°) 


n\ 
\'i 




(9) 


(10) 

(±U) 




(~\ CM 

\ LZ ) 


(16) 


PG0050+124 


0.02 


0.0 


0.81 


0.07 


7.92 


5.74 


0.71 


0.51 


0.05 


1.05 


0.005 


0.09 


PG0838+770 


0.78 


103.9 


16.32 


7.39 


0.74 


0.50 


0.34 


0.22 


0.01 


0.11 


0.006 


0.05 


PG0844+349 


0.00 


0.0 


128.93 


12.75 


6.15 


4.12 


0.61 


0.41 


0.07 


1.29 


0.007 


0.13 


PG0953+414 


0.00 


0.0 


25.62 


41.88 


3.25 


2.17 


5.32 


3.55 


0.05 


0.86 


0.08 


1.41 


PG1001+054 


3.03 


155.2 


1.39 


0.99 


0.08 


0.05 


0.06 


0.04 


0.01 


0.07 


0.006 


0.05 


PG1004+130 


0.56 


12.0 


1.62 


2.81 


0.19 


0.13 


0.33 


0.22 


0.01 


0.17 


0.02 


0.30 


PG1116+215 


0.00 


0.1 


43.18 


37.38 


3.80 


2.54 


3.29 


2.20 


0.05 


0.87 


0.04 


0.75 


PG1126-041 


4.16 


1.3 


3.71 


0.32 


1.40 


0.94 


0.12 


0.08 


0.03 


0.49 


0.002 


0.04 


PG1211+143 


0.27 


137.9 


179.59 


29.19 


3.21 


2.15 


0.52 


0.35 


0.15 


1.08 


0.02 


0.18 


PG1229+204 


0.00 


0.0 


1.81 


1.73 


3.20 


2.14 


0.31 


0.21 


0.06 


1.11 


0.005 


0.11 


PG1244+026 


0.00 


0.0 


156.25 


8.48 


3.68 


2.46 


0.20 


0.13 


0.02 


0.35 


0.001 


0.02 


PG1307+085 


0.60 


0.0 


30.92 


20.21 


1.46 


0.98 


0.96 


0.64 


0.05 


0.93 


0.03 


0.61 


PG1309+355 


0.52 


78.3 


12.54 


11.94 


0.65 


0.44 


0.62 


0.42 


0.03 


0.28 


0.03 


0.27 


PG1351+640 


0.00 


88.3 


3.87 


0.75 


0.56 


0.38 


0.11 


0.07 


0.03 


0.26 


0.005 


0.05 


PG1411+442 


610.01 


1.8 


4.29 


0.87 


1.83 


1.23 


0.37 


0.25 


0.03 


0.44 


0.005 


0.09 


PG1426+015 


0.03 


413.5 


40.22 


7.40 


8.58 


5.74 


1.58 


1.06 


0.35 


1.34 


0.06 


0.25 


PG1435-067 


0.00 


171.2 


9.81 


4.09 


2.82 


1.89 


1.17 


0.79 


0.11 


0.76 


0.05 


0.32 


PG1440+356 


0.00 


552.6 


107.79 


16.61 


3.39 


2.27 


0.52 


0.35 


0.17 


0.53 


0.03 


0.08 


PG1448+273 


0.00 


0.0 


129.83 


13.28 


2.97 


1.99 


0.30 


0.20 


0.02 


0.30 


0.002 


0.03 


PG1501+106 


0.27 


130.5 


483.37 


14.54 


13.10 


8.76 


0.39 


0.26 


0.63 


4.65 


0.02 


0.14 


PG1613+658 


0.29 


4.7 


103.88 


45.59 


4.98 


3.01 


2.18 


1.46 


0.11 


1.79 


0.05 


0.79 


PG1626+554 


0.29 


0.0 


80.88 


37.89 


4.49 


3.01 


2.10 


1.41 


0.03 


0.58 


0.01 


0.27 


PG2130+099 


0.56 


103.5 


99.01 


9.47 


3.28 


2.20 


0.31 


0.21 


0.17 


1.42 


0.02 


0.14 


B2 2201+31A 


0.25 


0.2 


90.30 


250.63 


4.07 


2.72 


11.30 


7.55 


0.12 


1.96 


0.33 


5.44 


PG2214+139 


3.80 


9.2 


10.07 


1.06 


3.25 


2.17 


0.34 


0.23 


0.10 


1.59 


0.01 


0.17 



Note. — Col.(l): Galaxy name. Col. (2): Intrinsic column density in units of 10 22 cm -2 . Col. (3): Ionization parameter from 
the PEXRIV component in erg cm s -1 . Col. (4): Absorption-corrected 0.5-2 keV flux for the blackbody component in units of 
10~ 14 erg s _1 cm~ 2 . The 2—10 keV flux of this component is <~10 3 smaller, and so is insignificant . Col. (5): 0.5-2 keV luminosity 
in units of 10 42 erg s _1 for the blackbody component. Col.(6)-(7): Absorption-corrected 0.5-2 keV and 2-10 keV flux for the 
power law component in units of 10 -12 erg s~ 1 cm -2 . Col. (8)— (9): Absorption corrected 0.5—2 keV and 2—10 keV luminosity of 
the power law component in units of 10 44 erg s _1 . Col.(10)-(11): Absorption-corrected 0.5—2 keV and 2—10 keV flux in units 
of 10~ 12 erg s _1 cm~ 2 for the reflected component. Col.(12)-(13): Absorption-corrected 0.5-2 keV and 2-10 keV luminosity of 
the reflected component in units of 10 44 erg s -1 . 
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Table 11. Quantifying the Errors in the Hardness Ratio Method 



Exposure 
(1) 



Hard 
(2) 



Soft 
(3) 



HR 

(4) 



r 

(5) 



(6) 



-Pb.5-2 keV 

(7) 



% crror . 5 _ 2 keV F2-10 keV % crror 2 _i keV 
(8) 



(9) 



(10) 



Input 








1.8 


1.77 




2.95 




100 


139.6±11.7 


453.2±21.2 


-0.53±0.05 


1.80t^ 


1.77 


0.0 


2.94 


-0.3 


75 


104.7±10.6 


339.3±18.3 


-0.53±0.05 


i-80±g;i| 


1.77 


0.0 


2.94 


-0.3 


50 


70.1±8.1 


225.6±15.3 


-0.53±0.07 


i-80±g;lS 


1.76 


-0.6 


2.93 


-0.7 


25 


35.0±5.9 


113.3±10.6 


-0.53±0.09 




1.77 


0.0 


2.94 


-0.3 


20 


27.9±5.1 


90.1±9.5 


-0.53±0.10 


1.80t°;» 


1.76 


-0.6 


2.92 


-1.0 


15 


21.0±4.6 


67.8±8.3 


-0.53±0.12 


i.8ot° : ? 7 


1.77 


0.0 


2.93 


-0.7 


10 


14.0±3.8 


45.2±6.6 


-0.53±0.15 


1 sn+ - 41 

i - SU -0.32 

i.8ot8;25 


1.77 


0.0 


2.93 


-0.7 


5 


6.9±2.8 


22.7±4.9 


-0.53±0.21 


1.77 


0.0 


2.94 


-0.3 



1.0 x lO^ 5 " 
5.0 x 10 20 
1.0 x 10 21 
5.0 x 10 21 
1.0 x 10 22 
5.0 x 10 22 
1.0 x 10 23 
5.0 x 10 23 
1.0 x 10 24 



15 
15 
15 
15 
15 
15 
15 
15 
15 



21.0±4.6 
21.1±4.5 
20.8±4.7 
19.9±4.4 
18.5±4.3 
12.2±3.6 
8.3±3.0 
1.7±1.3 
0.5±0.7 



66.4±7.8 
61.5±7.7 
55.6±7.6 
30.4±5.5 
18.2±4.4 
1.68±1.28 
0.2±0.45 
0.01±0.06 
0.001±0.03 



-0.52±0.12 
-0.49±0.12 
-0.45±0.13 
-0.21±0.14 
0.01±0.17 
0.76±0.34 
0.94±0.49 
1.00±1.13 
1.00±1.93 



4 - '8-0.26 
1 70+ ' 31 

1 62+ - 31 
,+0.27 
0.26 
+0.32 
-0.28 
-0.99 



1.13] 
0.71" 
-1.04 



1.73 
1.59 
1.43 
0.79 
0.49 
0.05 



-2.3 
-10.2 
-19.2 
-55.4 
-72.3 
-97.2 



2.96 
3.08 
3.12 
3.60 
4.24 
6.30 



0.3 
4.4 
5.8 
22.0 
43.7 
114.6 



Note. — Col.(l): Simulated exposure time in ksec. Col. (2): Average counts in the hard band from 1000 simulations. Col. (3): 
Average counts in the soft band from 1000 simulations. Col. (4): The hardness ratio. Col. (5): Spectral index derived from the 
hardness ratio method. Col. (6): Input column density from within the source in units of cm -2 . Col. (7): Observed 0.5—2 keV flux 
in units of 10~ 14 erg s~ 1 cm~ 2 assuming only Galactic absorption. Col. (8): The percent error of the simulated 0.5-2 keV flux 
relative to that of the input model. Col. (9): Same as (7), but for the 2-10 keV band. Col. (10): Same as (8), but for the 2-10 keV 
band. 
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Fig. 1. — The spectra of the PG QSOs with the best-fit power-law model. 
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Fig. I. — cont. 
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Fig. I. — cont. 
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Fig. 2. — Comparison of the reduced \ 2 values from the power-law models with those of 
P05 (top) and C06 (bo ttom) for all PG QSOs. In the left panels, the symbols represent 
the iNetzer et al.l (120071 ) FIR SED classification of each source where black circles represent 
strong FIR emitters, dark gray circles weak FIR emitters, light gray circles undetected FIR 
emitters, and open circles sources with unknown SED classifications. In the right panels, 
the symbols represent radio loudness. With the exception of two sources (PG 1211+143 
and PG 1501 + 106), the statistics suggest that the power-law models are as good, or better, 
descriptions of the quasar spectra as the P06 and C06 models. 
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Fig. 3. — Comparison of the spectral index for the power-law models with those from the 
C06 (left panels) and P05 (right panels) models. The meaning of the symbols is the same as 
that in Figure |2] for the SED classifications. The dotted lines are one-to-one ratios to help 
guide the eye. The indices based on the power-law and the P05 models tend to be steeper 
than those from the C06 model. Almost all of the indices are softer than the canonical value 
of 1.8 for AGNs but are still within the upper end of the range observed in other AGNs. 
Because the P05 model is not universal for the PG QSOs, the P05 figure is missing data 
from sources in their sample that do not overlap with our sample. 
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Fig. 4. — Comparison of 2-10 keV luminosity and the IR luminosity for the PG QSOs and 
the U/LIRGs. The horizontal axis is the IR luminosity and the vertical axis is the absorption- 
corrected 2-10 keV luminosity divided by the IR luminosity. The symbols for the PG QSOs 
are the same as those in Figure [2] for the SED and the values are from the global modeling in 
§ 14.21 For the U/LIRGs, the open stars represent the values derived from spectral fitting and 
the skeletal stars represent the values derived from the HR method. The colors of the stars 
symbolize the effective optical depth (r e //) taken from V09a where red represents objects 
with the highest r e //, green intermediate r e //, and blue the lowest r e //. The black stars are 
values from the RBGS archive sample (median Ljr = 10 1L98 L Q ), where r e ff is unavailable. 
The magenta squ are is the absorption-corrected value for Mrk 273 derived from Suzaku data 
(ITeng et al.ll2009l ) linked with the value derived from only the Chandra/ XMM- Newton data 
to demonstrate the effects of improved absorption correction. The solid line is the average 
log (L 2 -iokev/ Lir) for the PG QSOs (~ -1.6) and the dotted lines represent 1-a. The log 
hard X-ray-to-IR luminosity ratios of the U/LIRGs span about three orders of magnitude. 



-52 - 




Fig. 5. — Comparison of the bolometric luminosity from V09a with the absorption-corrected 
2-10 keV luminosity. The symbols for the PG QSOs are the same as those in Figure[2]and the 
symbols for the U/LIRGs are the same as those in Figure IU In the left figure, the horizontal 
axis is the total bolometric luminosity. For the PG QSOs, the 2-10 keV to bolometric 
luminosity ratio is nearly constant across the full range of bolometric luminosities. The solid 
line is the average 2-10 keV to bolometric luminosity ratio for the PG QSOs and the dotted 
lines represent one standard deviation of the mean. The total 0.5-10 keV luminosity of the 
PG QSOs are ~ 0.5 — 11% of the bolometric luminosity. Nearly all of the U/LIRG values 
fall below this trend. The right panel plots the same 2-10 keV luminosity as a function of 
the AGN bolometric luminosity (Xboi.AGN = /AGN-^boi) from V09a. The lines are the same 
as those in the left panel. Once again, a tight range is seen among the PG QSOs and most 
U/LIRGs fall below the relation. 
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bol. ACN/ ^Edd 



Fig. 6. — We compare the X-ray and Spitzer-dehved Eddington ratios for the U/LIRGs 
(same symbols as Figure HJ) and the PG QSOs (same symbols as Fi gure El) . The b l ack hol e 
masses used to calculate the Eddington luminosity are taken from IVeilleux et al.l ( 12009bl ). 
The PG QSO values determined by the two methods are linearly related. The dotted line 
represents this linear relationship where log ERx = (0.97 ± 0.52) log ERm — (1.88 ± 0.51) 
with R 2 = 0.43. This correlation is significant at the 99.97% confidence level. Unlike the 
quasars, most U/LIRG values do not follow this relation. 
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Fig. 7. — Comparison of the goodness of fit between the different models for the PG QSOs 
as a class. The horizontal axis is energy in the observer's frame and the vertical axis is 
the residuals in terms of sigmas. These are residuals for each spectrum plotted in the same 
frame to show, generally, where each model fails to describe the quasar spectra. Four different 
models are compared for the continuum of the X-ray spectrum. The residuals are binned 
to at least 10-sigmas for display only. Clockwise from the top left, the models presented 
are (1) MEKAL plus two partial covering absorbers (PCF; xl ~ 1-47), (2) a redshifted 
blackbody plus a neutral reflection model (PEXRAV; xl ~ 1-29), (3) a redshifted blackbody 
plus an ionized reflection model (PEXRIV; xl ~ 1.18), and (4) a redshifted blackbody plus 
a more complex and blurred ionized reflection model (REFLION; xl ~ 1-38). All of the 
above models include absorption by the Galaxy and an underlying power-law model for the 
AGN component. The PEXRIV model is the preferred model based on the fitting statistics 
and the residuals. Note that there is an absorption feature between 0.5 and 0.7 keV in the 
observer's frame. For the redshift range of our objects, this is consistent with the atomic 
transition of OVII and OVIII (E ~ 0.7 keV). 
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Fig. 8. — This is the same figure as the PEXRIV panel in the previous figure except with 
an additional redshifted absorption edge at 0.681qoo3 keV with r ~ 0.311qq3 (xt ~ 1-13). 
The absorption feature is consistent with the atomic transitions of O VII or O VIII. The 
addition of the absorption edge component improves the fit significantly, with A% 2 ~ 338 
for a change in two degrees of freedom. The large absorption feature still seen near 0.7 keV 
in the residuals is due to one source, PG 0050+124. 
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Fig. 9. — Plots of the absorption-corrected blackbody (or, equivalently, the soft excess) 
luminosity versus power-law luminosity from the best-fit PEXRIV model to the PG QSOs. 
The meaning of the symbols is the same as that in Figure [2j The solid line is the linear 
regression fit to the data except for PG 0050+124, the outlier near (44, 41). This source 
has an unusually strong absorption edge that the global model does not properly account 
for (see Figure [S]). The correlation is significant (>99.99%) with R 2 = 0.65 for logL bbody = 
(1.04 ± 0.28) logLp L - (2.69 ± 12.28). The slope of the line implies that L PL oc L bbody . This 
result suggests a direct physical link between the soft excess emission and the power-law 
component. 
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Fig. 10. — A histogram of the th ermal tempera t ure of t he accretion disk as derived from the 
photometric black hole mass of IVeilleux et al.l ( I2009bl ) and the 2-10 keV luminosity. The 
scatter is small and the histogram peaks near 1 eV with a standard deviation of 0.4 eV. 
These are only estimates of the temperature because the 2-10 keV to bolo metric correction 
for AGNs is uncertain. The predicted disk temperatures for the quasars are too low to 
explain the observed soft excess temperatures. 
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Fig. 11. — AGN contribution to the bolometric luminosity of the QUEST ULIRGs (same 
symbols as Figure Hj) and PG QSOs (same symbols as Figure [2]) derived by V09a against 
the logarithm of the 2-10 keV to bolometric luminosity ratio. The AGN% is the average 
of the six different methods (e.g., fine structure line ratios, PAH equivalent widths, mid- 
infrared continuum colors, and mid-infrared to far-infrared flux ratios) of measuring the 
AGN contribution by V09a. The linear relationship (dotted line) is significant (> 99.99%) 
and/? 2 = 0.64: AGN% = (21.11±3.45) log(L a _i fc e y/I'w) + (125.89±11.04). The correlation 
between the Spitzer-derived AGN% and the hard X-ray to bolometric luminosity ratio implies 
greater 2-10 keV flux with higher AGN contribution as one would expect. The Spitzer- 
derived AGN% are likely more uncertain in objects with intermediate and high T e ff, but 
this does not seem to be the origin of the scatter in the relation. Instead we favor intrinsic 
variation in the 2-10 keV to bolometric luminosity ratio for pure AGNs or unsuspected 
obscuration of the hard X-rays. See discussion in § 15.21 for more detail. 




Fig. 12. — A comparison of the 2-10 keV to bolometric luminosity ratio (our proxy for 
fractional AGN contribution to the bolometric luminosity; left) and the X-ray determined 
log (L/2-iokeV I Lem) (right) with some of the key physical properties of the U/LIRGs and 
PG QSOs: 25-to-60 fi dust temperature (top), optical spectral type (middle), and interaction 
class (bottom). The symbols for the PG QSOs are the same as those in Figure |2] for the 
SED and the symbols for the U/LIRGs are the same as those in Figure HI In the bottom two 
rows, the line connects the median values for each type/class of objects. The infrared-warmer 
objects have distinctly higher hard X-ray to bolometric luminosity ratios. The more Seyfert- 
like ULIRGs and the more advanced mergers (IVb and V) also tend to have a stronger AGN 
component, and the PG QSOs extend these trends. These trends are similar to those found 
in the analysis of the Spitzer data by V09a. 
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Fig. 13. — The starburst bolometric luminosity is plotted as a function of merger stage. 
O n the horizontal axi s, we have binary U/LIRGs (interaction classes Illab-IVa, phase D 
of lHopkins et al.ll2008l ). single-nucleus U/LIRGs (interaction classes IVb-V, phase D), FIR- 
strong PG QSOs where star formation is still present but the black hole dominates the 
feedback process (the "blowout" stage, phase E), and the FIR- weak PG QSOs where star 
formation has stopped (phase F). The open circles trace the mean total bolometric luminosity 
of the objects from both the QUEST and RBGS samples, the filled circles show the mean 
Spitzer-derived starburst bolometric luminosity from V09a for the QUEST sample only, and 
the filled squares represent the mean absorption-corrected 2-10 keV luminosity for both the 
QUEST and RBGS samples multiplied by a factor of 50 to be placed on the same scale as the 
other quantities. The error bars represent one standard deviation in each category of objects. 
This figure suggests that: (1) the contribution of the starburst to the total bolometric 
luminosity decreases as the merger progresses, (2) there is essentially no difference in AGN 
power between the FIR-strong and FIR-weak PG QSOs, (3) the growth of the AGN occurs 
most rapidly after coalescence, and (4) the 2-10 keV to bolometric luminosity correction for 
these AGNs is ~50, the normalization factor between the infrared and X-ray values. The 
large error bars on the X-ray data may be attributed to poor absorption correction or a 
broad intrinsic distribution of the hard X-ray to bolometric ratios among pure AGNs. 
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Fig. 14. — The fractional errors of the HR method from the average of 1000 simulations. The 
output photon index (black circles), 0.5-2 keV flux (red squares) and 2-10 keV flux (blue 
triangles), both corrected for Galactic absorption from the simulated Chandra spectra, are 
plotted against the input intrinsic column density. It is clear from the plot that the photon 
index derived from a single power law model deviates from the input value of 1.8 when the 
source column density is > 10 21 cm -2 , becoming flatter (Thr ~ 0.71 at N# of 1 x 10 22 cm -2 ). 
The 0.5-2 keV flux, more readily affected by absorption than the 2-10 keV flux, follows the 
same trend as the photon index. On the other hand, the hard-band flux is stable up to 
Nh ~5 x 10 21 cm~ 2 . These results of the simulations demonstrate that a flat spectrum is 
an indication of obscuration in the source. 
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Fig. 15. — A comparison of the values of T derived from the traditional spectral fitting 
method (horizontal axis) and the HR method (vertical axis). Included are the 13 objects 
from the QUEST sample (circles) that have enough counts for spectral fitting and also have 
been determined to have a power-law component in their spectra. We add to these 13 more 
nuclei from the RBGS sample (squares). The error bars are at 90% confidence level. The 
dotted line is the line of equality to help guide the eye. For most of the objects (filled 
symbols), the values of T derived from both methods are consistent with each other to 
within the errors. For the six sources with open symbols, the hardness ratio method severely 
underestimates T. All six objects have Njj > 10 22 cm~ 2 . This plot demonstrates that the 
hardness ratio method is a good estimator of the spectral properties of these faint sources as 
long as the column densities are < 10 22 cm -2 , consistent with the results of our simulations 
(Table H]). 
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Fig. 16. — A comparison of the 0.5-2 keV (left) and 2-10 keV (right) flux values of the 
U/LIRGs derived from spectral fitting (horizontal axis) and the HR method (vertical axis) 
for objects with enough counts for spectral fitting using c-stat. The flux values are in units 
of 10 -14 erg s -1 cm -2 . The bottom figures are close-up views of the boxed regions in the top 
panels. The solid line is a line of equality with each dotted line representing a 10% deviation 
from the spectral fitting values. Most of the HR values are within 50% (the shaded regions) 
from the fitted values. The HR method tends to overestimate the fluxes, especially at 2- 
10 keV, when obscuration is high. The median values for FHn/Ffi ts are ~1.2 and 1.3 for 
the soft and hard bands, respectively. This is an indication that many of these objects are 
obscured (see Table ITTT) . 



